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Abstract 



We consider systems of non-relativistic, interacting electrons at finite 
density and zero temperature in d = 2, 3, ... dimensions. Our main concern 
is to characterize those systems that, under the renormalization flow, are 
driven away from the Landau Fermi liquid (LFL) renormalization group 
fixed point. We are especially interested in understanding under what 
circumstances such a system is a Marginal Fermi liquid (MFL) when the 
dimension of space is d > 2. The interacting electron system is analyzed by 
combining renormalization group (RG) methods with so called "Luther- 
Haldane" bosonization techniques. T RG calculations are organized as a 
double expansion in the inverse scale parameter, A -1 , which is proportional 
to the width of the effective momentum space around the Fermi surface 
and in the running coupling constant, g\, which measures the strength of 
electron interactions at energy scales ~ Vp £ F . For systems with a strictly 
convex Fermi surface, superconductivity is the only symmetry breaking 
instability. Excluding such an instability, the system can be analyzed by 
means of bosonization. The RG and the underlying perturbation expan- 
sion in powers of A -1 serve to characterize the approximations involved 
by bosonizing the system. We argue that systems with short-range inter- 
actions flow to the LFL fixed point. Within the approximations involved 
by bosonization, the same holds for systems with long-range, longitudinal, 
density-density interactions. For electron systems interacting via long- 
range, transverse, current-current interactions a deviation from LFL be- 
haviour is possible : if the exponent a parametrizing the singularity of the 
interaction potential in momentum space by E(|p|) ~ ii is greater than 
or equal to d — 1, the results of the bosonization calculation are consistent 
with a MFL. 



1 Introduction 



In this paper, we consider systems of non-relativistic electrons at finite density 
and zero temperature in d = 2,3, ... dimensions. The interactions between elec- 
trons are described by two-body potentials or by current- current interactions. 
The Cooper channel which drives the BCS instability is turned off; (e.g. by 
assuming that the Fermi sphere of the non-interacting system has a suitable ge- 
ometry |, |). 

Our main concern in this paper is to characterize those systems that, under 
the renormalization flow, are driven away from the Landau Fermi liquid (LFL) 
renormalization group (RG) fixed point. More concretely, we are interested in un- 
derstanding under what circumstances such a system is a Marginal Fermi liquid 
(MFL) when the dimension of space is d > 2. This problem comes up, for exam- 
ple, in the study of single-layer quantum Hall fluids at filling fractions v = ^, |, ... 
and, perhaps, in the theoretical description of materials related to anisotropic 
HT C superconductors (see || f|, |5j, |6| and references therein) . 
A one-dimensional MFL (or Luttinger liquid) at zero temperature can be char- 
acterized by the property that the electron propag ator falls off like \x\-( 1+, i\ at 
large distances \x\, for an exponent r\ > that depends on the electron-electron 
interaction and characterizes a new RG fixed point. When rj vanishes the system 
is a LFL . In the limit of large distance scales and low frequencies (scaling limit) 
the properties of a LFL are identical to those of a free system of non-interacting 
electrons (up to a renormalization of the residue of the one-particle pole and of 
the Fermi velocity). In one dimension, MFL correspond to a line of RG fixed 
points in the space of effective Hamiltonians (or actions) containing the fixed 
point corresponding to the free system. 

In d > 2 dimensions, we define a MFL to be a Fermi liquid with an electron 
propagator falling off more rapidly than the free electron propagator by at least 
a fractional inverse power of the distance between the arguments, but not ex- 
ponentially in the distance. Contrary to the superconducting instability, the 
instabilities leading to a MFL are not accompanied by symmetry breaking, and 
there is no energy gap in the excitation spectrum of such systems. Yet, Landau's 
picture of non-interacting quasi-particles does not apply to the physics of Lut- 
tinger liquids. 

Historically, MFL were discovered in the context of one- dimensional systems of 
interacting electrons , where they appear naturally, for a large class of two-body 
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interactions. The experimental observation that the normal phase of anisotropic 
HT C superconductors exhibits many non-conventional features, incompatible with 
LFL theory, leads to the question whether electron-electron interactions can drive 
a system in two or more dimensions away from LFL's. The same question arises 
in the context of single-layer QH fluids at filling fractions v with even denomina- 
tors. 

Recently, there have been many investigations of this problem. Besides numeri- 
cal studies, two analytic techniques proved useful in attempts to understand it: 
the RG method, involving an expansion in an inverse scale parameter A -1 ||, 
which has led to a wealth of rigorous results |], || ||, |i~0f , and the bosonization 



techniques [|TTJ proposed, in this context, by Haldane fl"2| . 

In this paper, we study the stability and instability of the LFL by using the 
same two methods. We propose to combine them to clarify the picture of Fermi 
liquids. We implement the bosonization technique in a manner elucidating its 
standing in terms of the fermionic perturbation expansion and rendering the cal- 
culation of the electron propagator at large distance scales quite transparent. 
Our techniques apply to a broad class of two-body interactions. (We presented 
a preliminary account of our results at the 1994 Les Houches summer school 



"Fluctuating Geometries in Statistical Mechanics and Field Theories", fll7| ) 
Next, we summarize the contents of the various sections of this paper. 



In Sect. 2, we review the RG method for non-relativistic electron systems at 
finite density and zero temperature in a d+ 1-dimensional, euclidean space-time, 
with d > 1. The underlying perturbation theory is organized as a double ex- 
pansion in an inverse scale parameter A -1 and in the (dimensionless) running 
coupling constant, g\, of the two-body interactions. 

The inverse scale parameter A -1 is proportional to the width of a shell Q\ around 
the Fermi surface. Having integrated out the electron modes with momenta lying 
outside the shell fl\, the effective action, S\ , for the remaining modes describes 
the physical properties of the electron system at distance scales ~ or energy 
scales ~ Vp ^ F , kp and vf being the Fermi momentum and velocity, resp. . The 
running coupling constant g\ measures the strength of the two-body interactions 
in S\. 

Given S\, one further reduces the effective momentum space around the Fermi 
surface to fly, with A' = MX , M > 1 , and calculates the corresponding effective 
action Sy. This calculation is organized as a perturbation expansion in the pa- 



3 



rameters g\ and A -1 , assuming that they are small. Rescaling the resulting system 
by a factor M, we obtain a system that is similar to the initial one, and we can 
compare the relevant and marginal parameters characterizing the corresponding 
effective actions. The RG method consists in executing these transformations 
iteratively, with the aim of deriving the scaling limit of the system, as A — > oo. 
The underlying calculations can be interpreted in a transparent way by decom- 
posing the shell Q\ into iV ~ A d_1 boxes B Ui {\) , i — 1, N , with sides of length 
~ t. The d-dimensional unit vectors Ui , i = 1, N , point to the centers of the 
corresponding boxes. This decomposition implies a decomposition of the electron 
field ^ (where \&" stands for \&* or \P) into N components ip^. , % — 1, ...,iV, 
whose Fourier modes have support in the corresponding boxes. A component 
field describes quasi 1 + 1-dimensional, "relativistic" electrons moving along 
the direction Ui with velocity Vp. The decomposition of the electron field into 
such component fields reproduces the electron propagator to leading order in an 
expansion in j-. In the non-interacting system, the subsystems along the different 
directions are independent. Two-body interactions couple the subsystems. After 
rescaling, each interaction process is suppressed by a factor Hence, the inter- 



acting system can be analysed by using " large- iV" expansion techniques, || |9], |14 
For systems with a strictly convex Fermi surface, superconductivity is shown to 
be the only symmetry-breaking instability that can develop in the system. We 
propose to study electron systems which do not undergo such a symmetry break- 
down - i.e. whose Cooper channel is turned off - by means of "bosonization" . 
The framework of the RG will serve to characterize the approximations involved 
by "bosonizing" the system. 

In Sect. 3, we introduce the bosonization technique by calculating the scaling 
limit of the effective gauge field action (the generating function of connected cur- 
rent Green functions) for a non-interacting electron system. The effective gauge 
field action is obtained by coupling an external gauge field to the electron system 
and by integrating out the fermionic degrees of freedom. 

The calculation of the scaling limit of the effective gauge field action is reduced 
to the calculation of the gauge field action of a family of independent Schwinger 
models, because the d- dimensional, non-interacting electron system decomposes 
- to leading order in an expansion in the inverse scale parameter A -1 - into inde- 
pendent subsystems of quasi 1 + 1 dimensional, "relativistic" fermions, one along 
each direction [u] = {uj, —uj}. 

Gauge invariance must hold for each subsystem - i.e., in each direction - sepa- 
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rately, and implies local conservation laws for the associated quasi 1 + 1 dimen- 
sional current densities , A = 0,1. These conservation laws imply that each 
current density jfo can be expressed as a derivative of a bosonic field ip^y It turns 
out that, for the non- interacting system, these fields are massless and gaussian, 
propagating only along the direction [uj\. 

By taking the effective action S\ determined in the context of an RG analysis as 
an input, we can apply the bosonization technique to interacting systems, too. 
For systems whose Cooper channel is turned off, the electron-electron interac- 
tions are described - to leading order in the inverse scale parameter A -1 - by an 
expression quadratic in the currents jfa . By replacing the fermionic currents jfa 
by the corresponding bosonized expressions, we obtain a gaussian bosonic theory. 
Because the theory is gaussian, it can be studied more easily than the original 
fermionic theory where the interaction term is quartic in the electron fields. 
In the calculation of the effective gauge field action, this technique is shown to 
reproduce the leading order of a (fermionic) perturbation expansion in powers of 
A -1 . The leading contributions coincide with the ones of an RPA approximation. 
We determine the scaling limit of the effective gauge field action for systems with 
short range interactions and discuss the extension of the bosonization method to 
systems with long range interactions in three examples : systems with longitudi- 
nal density- density or transversal current-current interactions, and " Tomographic 
Luttinger Liquids" . 

In Sect. 4, we apply the bosonization technique to the calculation of the elec- 
tron propagator at large distance and time scales. We start from an effective 
action S\ at energy scales ~ and the corresponding decomposition of the 
electron fields ^ into N quasi 1 + 1 dimensional components ipl,. We replace 
the action S\ by its bosonized version, introduced in Sect. 3, and , for each ray 
[uj] = {uj, —lu} , we express the pair i/j^i/j^ in terms of a bose field ip^}. This 
is accomplished by applying the well-known bosonization formalism for 1 + 1 di- 
mensional relativistic fermions . However, one has to cope with a subtlety arising 
from the dependence of the quasi 1 + 1 dimensional electron fields ip^ on the 
components of the momentum perpendicular to the direction [uj] . Again, the 
approximations involved can be characterized in the context of a (formal) per- 
turbation expansion in powers of the inverse scale parameter A -1 . Special care is 
taken to discuss the implications of linearizing the Fermi surface inside the boxes 
5 W (A). 

We verify that the electron propagator of a system with short-range interactions 
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tends, for large arguments, to the standard Landau Liquid form. Because of 
screening - which is reproduced by bosonization - the same result holds for sys- 
tems with long-range density- density interactions. 

For systems of electrons interacting via long-range, transverse current-current in- 
teractions, we find the possibility for a deviation from LFL behaviour, depending 
on the exponent a which characterizes the singularity of the interaction poten- 
tial in momentum space V(p) ~ • The critical value for a is d — 1. For 
< a < d — 1 and A —>■ oo, the electron propagator has the standard LFL form, 
whereas, for d — 1 < a < 2 and A — > oo,we argue that it shows MFL behaviour. 
In the second case, in order to obtain a theory which is form-invariant under scale 
transformations, the parameters of the effective action S\, used as an input of 
bosonization, must be functions of the scale parameter A. A resulting consistency 
condition determines the expected flow of the parameters of S\ under RG trans- 
formations. 

Our analysis yields a view somewhat complementary to those arrived at in || |], 
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2 Effective Action on Large Scales and Fermionic 
Perturbation Theory 



The action of a system of non-interacting, non-relativist ic electrons in d+ 1 space- 
time dimensions is given by 



= [dtd ( 



x 



t) (ido + fj,)V(x, t) 



2m 



where m denotes the bare electron mass and the chemical potential, /x, specifies 
the mean electron density n° . We choose units such that h — 1. 
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In a functional integral quantization, the electrons are described by two indepen- 
dent, two- component Grassmann fields t) = ( ) an< ^ ^*{x,t) = 
( tp*(x,t) ip*(x,i) J . The arrows |, j, stand for "spin up" and "spin down". 

An ionic background and interactions between electrons are described by ad- 
ditional terms in the action of the form 



dtd d xU(x)^*(x)^(x) 



and 



S v = —J dtd d x J dsd d y m*(x)m(x)V{x-y) 5{t - s) V*(y)V(y) , 

resp., where U is a one-body potential and V a two-body potential. The total 
action is given by 

S = S° + S u + S v . (2) 

Expectation values of functionals, J- (iff*, iff) , of the iff and 1 J r * fields are calculated 
by evaluating Berezin- Grassmann integrals : 

( T(iff* iff) ) = 1 1 1 { 1 1 (3) 

\j-y*,* )/fl JV(iff*,iff) e^C**.*;/") ' KJ 

We work in the Grand Canonical Ensemble, where the chemical potential \i is 
held fixed. This implies that, in general, the mean electron density of the 
interacting system does not coincide with the mean electron density n° of the 
free system. 

In order to calculate expectation values of the form in eq.(|]), the total action S 
is split into a quadratic part, S 2 = S° + S u , and a quartic part S v . One expands 
the exponential e lS in a power series and calculates the expectation values of 
the resulting polynomials in iff* and with respect to the Gaussian Berezin 
integration determined by S 2 , by using Wick's theorem. 

In this paper we consider only systems invariant under translations and rotations 
of space, i.e., the background is described by a constant potential U(x) = U that 
can be absorbed in a redefinition of the chemical potential ("jellium model"). 

For technical convenience, we analyse the system in the euclidean region 
reached by analytic continuation in the time t to the halfplane 9m[t] > and 
setting x = it (Wick rotation). For a system at a finite temperature T, the 
Zo-variable in the euclidean action is integrated over the interval [— | , | ] , where 
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(3 is proportional to the inverse temperature T _1 , and anti-periodic boundary 
conditions are imposed at x = ±|. In this paper, however, we only consider 
systems at zero temperature, (f3 — > oo). 

Then the covariance of the Gaussian integration - i.e., the unperturbed electron 
propagator - is given by 

r r e -ipo(xo-yo)+ip(S-y) 

= 8^ Jjj* J K J*p m _ £ o m . (4) 



where = ^ — p. is the energy of a free electron, and the indices a, (3 label 

the spin orientations. We use the short-hand notation / dp(-) := / (■)• 

In the ground state of the unperturbed system, all one-particle states with 
wave vectors p satisfying \p\ < kp = \\J1m\i are occupied, and the electron 
density is given by 

nl = 2jd d p9(k F -\p\). (5) 

The Fermi surface , S%~ 1 , is defined as the surface of this sphere and the Fermi 
wave number, fcjr , sets the fundamental momentum scale of the system. 
For the main results and conclusions of this paper to hold, the spherical symmetry 
is not essential. Our analysis applies as long as the Fermi surface is strictly con- 
vex. The situation changes, however, radically for systems of electrons hopping 
on a square lattice, at half filling, where the Fermi surface contains flat, parallel 
faces giving rise to "nesting phenomena". In addition to the superconducting 
instability, other instabilities like charge- or spin- density-wave instabilities can 
occur, and their interplay can lead to rather complicated phenomena. Such sys- 
tems can be analyzed by using methods similar to those described in this section, 
but their properties are not yet fully understood. 

We are interested in universal large-scale and low-energy properties of electron 
systems, i.e., in the so called scaling limit of such systems. The scaling limit can be 
constructed by using renormalization group (RG) techniques which are based on 
successively integrating out the modes of the electron fields , \P corresponding 
to wave vectors far from the Fermi surface, with the aim of deriving an effective 
action for the modes close to the Fermi surface; (the energy vanishes for modes 
whose momenta lie on the Fermi surface) . 

An alternative way to obtain some "non-perturbative", large-scale and low-energy 
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information about an electron system consists in using the so called "Luther- 
Haldane" (LH) bosonization technique [|ll], [12|, |13[ . However, for electron systems 



in more than one space dimension, LH bosonization is not an exact method, i.e., it 
does not exactly resum the perturbation expansion in the coupling constant, g, of 
the quartic electron-electron interaction . The aim of this paper is to incorporate 
the bosonization technique into the systematic framework of the RG, in order to 
estimate the effects of the approximations involved by " bosonizing" an electron 
system in more than one space dimension. 

In general, RG calculations for interacting electron systems require assuming 
that the coupling constant g of the quartic electron-electron interaction is small. 
For one-dimensional systems, the calculations have been carried out a long time 
ago (for a review, see [0]). In more than one space dimension, the situation 
is more complicated. However, during the last few years, there has been sub- 
stantial progress accomplished by introducing the inverse scale parameter , A -1 , 
proportional to the width of the effective wave vector space around the Fermi 
surface, as a supplementary expansion parameter, ||. In this way, the dominant 
contributions to the scaling limit are obtained in a natural manner . 

Here we sketch a version of the RG method for electron systems in more 
than one space dimension based on a double expansion in A" 1 and the running 
coupling constant, g\, which measures the strength of electron interactions at 
energy scales ~ for original results, see [| |10[ , and for reviews, 

13 , 17 1 . (An alternative way would consist in artificially introducing S different 



species of electrons and to organize the perturbation expansion in powers of A 1 
and S^ 1 or S, for S — > oo or S — > 0, resp. ; see Hj.) 

We adopt a Wilson-type formulation of the RG and accomplish the mode 
reduction in momentum space. The electron fields \l/*(x), ty(x) are expressed in 
terms of their Fourier modes 

#(p) = J d d+l x e ipx V{x) **(p) = J d d+1 x e~ ipx ^*(x) , 

with px := p Xq — px . 

We assume that, for a sufficiently small value of the dimensionless coupling con- 
stant g, there exists a large scale factor Ao ^> 1 (with Ao — > oo, as g — > 0), such 
that, in a first step, the integration over the electron modes with 
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momenta p outside the shell 

p- tip- 
Pi 



Q Xo := {peR d , 



< 




of width ^ around the Fermi surface leads to an effective action that has essen- 
tially the same form as the original action S. More precisely, the effective action 
for the remaining modes with momenta inside the shell Q\ is given by 



Sj[(**(p),*(p)) = -log/ v(*%p')Mp))e^^ + aSV ^ 

+ cte. . (6) 

We assume that 

s£(**(p),*(p)) ~ s(**(p),i(p); 



In other words, for sufficiently small coupling constants (< ^) , possible insta- 
bilities develop only at energy scales smaller than Vp-^ in Q, it has been 
proven rigorously that this is true for a two-dimensional system with short-range 
interactions. 

Given S)~ , we wish to determine the effective action, 5^, on a lower energy 
scale ~ Vp^- , where Ai = MXq , M > 1 , by integrating out the modes in the 
shells Q\ \ Q\ 1 . Functional integration leads to a perturbation expansion for 
SJ^ with a fermion propagator determined by the quadratic part of . 

In order to illustrate characteristic features of the resulting perturbation the- 
ory, we determine the unperturbed electron propagator G°(x) , defined in eq.([|), 
for large values x = A£ of its arguments, where £ ~ and A —>■ oo . One finds 



F 



27T J Si 



e -i(poA£o-p||Atj£) 



r lpo - VFpil 



(7) 



cf. [0, where a; denotes a <i— dimensional unit vector (Co £ ), and the Fermi 
velocity vp = ^ is given by the linearization of the energy function e° (kpu+p) = 
vp\p ■ uj\ + 0(p 2 ) around the Fermi sphere . For large arguments (~ A) of the 
electron propagator, only momenta inside the shell are important. 
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One can subdivide the shell Q x into N = (^) d_1 Vol(^; 1 ) disjoint, congru- 
ent blocks B u) (i)(X) with sides of length ~ &f , centered at the points kptu^ , % = 
1, ...,N , on the Fermi surface, as indicated in Fig.l .The endpoints of the unit 
vectors , % = 1, N , form a regular grid on Sf~ l . 




Fig.l 



This discrete decomposition of the shell Q\ into iV boxes B u) (i)(X) yields a 
formula analoguous to eq.(|7]) for the large distance behaviour of the electron 
propagator : 



N 



8=1 



dpo /_ ■ 

R JBj t) (X) 

e -*(poA?o-pA^ 

ip - £°(A;f^ w +p) 



l+0(-) 



with S wW (A) := {q- k F <2® , qe B w(i) (X)} . 

The energy function e° of the unperturbed system can be written as 



(9) 



where the momenta p are supposed to lie in the box B u (i)(X) , with p\\ : = 
P|i a; , pi! = (uj ■ p) , and p± := p — f>\\ ■ Compared to the first, linear term, the 
quadratic contributions are of higher order in ~ , and one is tempted to neglect 
them. Neglecting the dependence on the perpendicular momenta, p±, amounts 
to replacing the piece of the Fermi surface, S^ 1 n B U} (i ) {X) , in the box B u) ( i )(X) 
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by a plane, cf. Fig 2 . 

This approximation is harmless as long as p\\ 3> (^f) 2 • I.e., given the decom- 
position of the shell fl\ into boxes B u (i) (A) , one can neglect the corrections to 
the linear part of the energy function e° in the propagator, eq.(|j), as long as one 
only integrates out modes with momenta py ^> (^) 2 . 




Fig.2 



After having linearized the energy function e° eq.(§) can be reproduced by 
decomposing the electron fields ^f* , \& into N independent components r/>* w , ip^i) 

N N 
i=l i=l 

(10) 

with propagators 

-i(fco(^ -J7o)-fc(C-^)) 
^Jfc- — ^- . 

(i) i£; — vpu^k 

Here, "~" stands for "equal to leading order in | " and B^i) '■= JB W (<)(1) . 
The Fourier modes of the component fields w (£) , i.e. (i) (k , ■ ) , with = 



<Vuo(0€(ofa)> 



R 



dk 
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\I/or\l/* , have support in the boxes B^w with sides of length ~ kp , and their 
propagators are given by 



y wW) uU ) W d+1 s (d+1) {k - k') =i 1 o (k) 



(11) 

where 

1* 0) := I I ' ^ GSwW 
B -w v ; \ , else 

Note that the momenta in eq.flTTD are related to the momenta p in eq.(§) by the 
scale transformation p — f . Whereas the electron fields ^(x = A£) are functions 
of the physical coordinates (x), the component fields ip (»)(£) are functions of the 
rescaled coordinates (£). In the rescaled system, the Fermi surface has the radius 
A/cf, and the boxes B u (i) cover a shell Q of thickness fcp around the Fermi surface. 
Eq.(|IIp shows that, in order to describe the large scale physics of the unperturbed 
system, the electron fields \l/ can be decomposed into N ~ A d_1 components 
VC> V'w which propagate only in the direction perpendicular to the Fermi surface. 

We now return to the interacting system. After having integrated out the 
modes with momenta in R d \ Q\ , we can apply the decomposition (|T0|) to the 



remaining electron modes in the effective action (p) , & (p)^J , the error 

being of order i. According to the assumptions stated above, the effective action 
Sq S = at an energy scale ~ then has the form 

Sf = + 5S* + S* + " higher-order terms " , (12) 



with 



r-2 
°0 



EE / ^yiWfe-wll^) (13) 
a, 0=1,1 Jl «> z o 

~ EE / ^ A ^ ^(fc) (14) 



and 



*o 4 - EE / / a d+1 k^ 



1 A;i_ 

1 r 

^0 W1,...,£J4 CT ,0" 



o^Cw^i) C 4 .(^ {4) )C3.'( fc(3) )^'(^ (2) )^ 1 .(fc (1) ), (15) 
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where u := {ui, d; 2 , u 3 , u; 4 } , k : = {^1,^2,^3,^4} , and I u stands for the inte- 
gration domain R x . 

For small values of #Aq, the parameters v F q, Z ,gQ' a (o> ; -^k) are renormal- 
ized only weakly with respect to the parameters v F , Z = 1, gV in the original 
action. The set of dimensionless coupling functions g^" ( u ; -±k) is related to 
the Fourier transform, gV(-j^), of the original interaction potential by 

k F - d gr\u; £4) « g ^ V (k F0 (u; 4 - cJi) + - . (16) 

The Fourier transform, V", of a short-range two-body potential V is smooth in 
momentum space, so that - to leading order in - we can neglect the dependence 



on the small momenta. Thus, the coupling functions g^" (u; ifc) in eq.(15) 
can be replaced by a set of coupling constants g^ a (uf) . 

For a long-range potential, whose Fourier transform V(p) is singular in momentum 
space, i.e., V(p) ~ t4^ , with a > , we set 

,i_d t^/ !^_J goS^V (k FQ (u A -uTi)) , fora} 4 ^a;i, 
' " ( - : ~ i <? ^ V (i (£ (4) - £ (1) )) , for tf 4 = ft • ( j 



We first restrict our analysis to short-range potentials; comments about long- 
range potentials will be made at the end of this section. 

In general, a quadratic term 5Sq of the form given in eq. (|T4T) is generated. 
It displaces the origin of the energy spectrum - i.e., the Fermi wave number - 
to \ k F0 m \ (k F — ^) . Under the condition that A ^- <C k F , we can 
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absorb this term in a change of the parallel momentum , k\\ — > ki = fey + Ao 
obtaining 

^ := S 2 + JSg = £ £ / ^ <a(*0 ( ik 'o- v ™uk>) ^ a {k>) . 

(18) 

2 

Thus the propagators obtained from are equal to the propagators determined 
by Sq , except for a small displacement of their support, B u — > . 

By "higher-order terms", in eq.(|12|), we mean contributions corresponding to 
higher orders in the Taylor expansion of the coefficient functions of the quadratic 
and quartic terms in the momentum variables, or contributions involving more 
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than 4 electron fields. Engineering scaling suggests that both types of contribu- 
tions are irrelevant, but we shall analyze these terms more carefully later in this 
section. 



As usual, one divides the terms in the action into relevant, marginal and 
irrelevant ones, in accordance with their scaling dimension. The exponent , v , 
characterizing the (leading) behaviour of an expression .F k) under 

scale transformations, i.e. , 

^(^A)(*).$(A)($);i) « A" ^ ($*(*), $(*:);*) , forA^oo, 

is called the scaling dimension of JF . The scaling dimension of the electron modes 
is | + 1, i.e., 

*Ui) = • (19) 



This scaling dimension is fixed by the requirement that the scaling dimension of 
the quadratic action S% is zero. 

Terms with scaling dimension greater than / equal to are called relevant / 
marginal, and terms with negative scaling dimension are called irrelevant, as 
they are suppressed by inverse powers of A , as A — > oo . 

The two terms of the quadratic action Sq , given by eq. (|T3D , are marginal. Con- 
tributions to Sq which arise from higher orders in the Taylor expansion of the 
energy function, e , in the momentum variables are not displayed, as they are 
irrelevant. The quadratic term 5Sq , cf. eq.(|l4"D , which causes a displacement of 
the Fermi wave number has scaling dimension 1 , i.e. , it is relevant. By holding 
the chemical potential /i fixed, the average electron density of the system - related 
to the Fermi wave number - changes during the RG iterations, and this requires 
a continual readaptation of the linearization point of the energy function. 
The quartic interaction , eq.(^) , has scaling dimension 1 — d , i.e. , it is ir- 
relevant in dimension d > 2 . However, there are of order N ~ Aq _1 different 



interaction terms, cf. eqs. (O) and (IH). It can happen that - for special ex- 



terior momentum configurations - the sum over the different interaction terms 
compensates the scaling factor -^rr (as discussed below ) . Moreover, engineering 
scaling arguments are only valid as long as the dimensionless running coupling 
constants remain small during the RG iterations. If an instability is developing, 
in the sense that some couplings diverge as A grows, our weak coupling analysis 
breaks down. The RG method will permit us to identify those processes that lead 
to instabilities. 
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Corrections due to higher orders in the Taylor expansion of the interaction poten- 
tial around the points kpoi^^ — ^ V) ) are sub-leading, cf. eq.flTBD . (Of course, 
this is only true for short-range interactions). 

Local terms in the action involving 2£ electron fields ( I > 2 ) have scaling dimen- 
sion — £d + (d + 1) ; they are irrelevant and can be neglected. However, non-local 
terms of the type shown in Fig. 3a (straight lines stand for electron propagators 
G u and wiggly lines for local 4- vertices) can appear which have scaling dimension 
[£ — 1)(1 — d) . These are tree-level contributions constructed exclusively out of 
£ — 1 local 4- vertices; (the non- locality arises from the inner propagator lines). 
Although they have a negative scaling dimension, they can contribute in leading 
order to the renormalization of some special, lower order vertices. This happens 
in cases where momentum conservation allows one to sum over the a;— labels of 
contracted pairs of incoming and outgoing propagator lines, leading to a contri- 
bution of 0(Ag _1 ) per contracted pair ( cf. Fig. 3b) . 



0>2 



UJ 2 



0J 2 



V 3 



^3 



Fig. 3a 




The influence of such higher order vertices to the RG flow of the parameters 



vp, Z, 5fj,, g(u) is analyzed in |16j ; they contribute in leading order in the inverse 
scale parameter A, but do not cause any qualitative change of the flow. For the 
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sake of simplicity, we shall ignore them in the following, because we assume that 
g is small. (However, in the RG teatment of systems with long-range interactions, 
they play an important role in obtaining the correct form of screening) . 



Given Sq , we attempt to calculate the effective action Sf at a larger scale 
Ai = MX > X . Assuming that the dimensionless quartic coupling constants 
go(uf) are small, i.e. |<7oGtOI < 9 ? we organise the perturbation theory in 
powers of g . The appearance of the small factor in front of the quartic 
interaction term, eq . (fTED , suggests to introduce the inverse scale parameter as 
a supplementary expansion parameter. This is also the parameter that controls 
the approximation of replacing the electron fields ^ by the component fields 
VCj^w; cf- ec f(0) • m order to obtain the effective action, £^ ff , on successively 
larger scales X n = M n X , n — > oo , we shall execute this procedure iteratively. 
This leads to the following formulation of the RG. 



In the iteration step n — > n + 1 , we lower the energy scale in the effective 
action S^ ff by a factor M which we choose to be an integer greater than 1 . The 
action 5^ ff is supposed to have the form specified in eqs. ( |12|) - fllED , with the scale 
parameter Ao replaced by A n = M u Xq , and the parameters Z q ,vfq, S/io, go(u) 
replaced by Z n ,vp n ,5^ n , g n (u_) . Then, the symbol stands for "equal to 
leading order in j-" . There are N n = (|^) d ~ 1 Vol( l S n ) boxes which cover 
the shell fl n of standard width k F around the spherical surface S n of radius 

The iteration can be continued as long as the running coupling constants g n {ui) 
remain small, i.e., |<7 n (^)| <3<Cl and the displacement, A n ^- , of the origin 
of the energy function from A n _iA;p n „ 1 to X n kp n := A n (&p n „i — ^ rL ) is small 
compared to ^j, cf. Fig. 4 . As in step 0, cf. eq. QT8D , this allows us to absorb the 
displacement of the chemical potential - described by 5S^ - in the term S 2 n by 
shifting the variables k\\ — > k'n = fcy + An^ 21 . This yields a displacement of the 
Fermi wave number from X n kFn-i to X n k Fn (« X n k Fn _i) . 
We denote the shell of width ^ around the spherical surface S n of radius X n kp n 
by Q n . In the step n — > n + 1, we eliminate modes with momenta k lying in 
Q n \Q n . This yields an effective action 



sf +1 ({A(k)}) « 

w -log / _ P{^(fc')ie-[ s " ff({ ^ (fc ' )}) ] + cte. , (20) 
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where {^} denotes the set of component fields with tp^ = -0* or t/j^ . In order 
to compare Sf +1 to Sf , we subdivide each box n Q n into M congruent 
boxes, B$ t , £—1, M , with sides of length ~ ^ (cf. Fig 4) . 




Fig. 4 



The resulting number of boxes is iV n+1 = M d ~ 1 N n . This requires a refinement of 
the decomposition (jT^) of the electron field. Note that, in this way, we implicitly 
take into account the curvature of the Fermi surface. Using eq.(^), we then 
rescale the momentum variables, k — > kl = Mk , with k' E B^ 1 ^ , where the 
boxes B^ +1 ^ have standard size. 

We end up with an effective action S^, 1 of the same form as S^ ff . The new 
effective parameters Z n+ i, vfu+i, 3fi n +i, g n +i{dl) can be expressed as functions 
of the set of parameters V n := {S^ n , Vp n , g n (u), A n } : 

Sfin+l 



A 



n+l" 



J n+1 



^-^(k;V n ) 

Zlr, 



lfc=0 



VFn+l 



1 

1 



1+i 







J n+1 



dk 

d 



lfc=0 



(21) 



VFn + -KT T>u{k;T> n ) 

okw 



\k=Q 
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and 

4^ - irfo&noi^ . (22) 

The functionals and T, which appear in these flow equations, turn out not to 
depend on the iteration step (to leading order in -£-). 

The functional is the self-energy and is obtained from the amputated one- 
particle irreducible (1PI) connected graphs renormalizing the propagator line 
Gu, n > an d r(w), the 4-vertex function, contains all 1PI connected graphs that 
renormalize the dimensionless coupling constant g n (ui) . We restrict our atten- 
tion to determining the functionals and r(u) to leading order in a double 




Fig.5 



In Fig.5 , we display diagrams contributing to the self energy S w in leading 
order in y- ("cactus diagrams") . To leading order in g , only the first one is 
retained. Straight lines stand for electron propagators G u>n given by 



and dots stand for the following combinations of interaction vertices —k F 
represented by wiggly lines : 



1— dgn{o£) 



X 



d-l j 




and 
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where momentum conservation requires that U\ + u 2 — <^3 + ^4 • 
In Fig. 5, a summation, YX') > over ^ £ is associated with every electron 

loop and compensates the factor — per interaction vertex ( for graphical rea- 
sons, the spin indices are omitted in the loops; we sum over all compatible spin 
configurations). One verifies easily that - for the case of short-range interactions 
- the contribution of the first diagram in Fig. 5 is a constant of order g . Using 



eq.fl2lD, this implies that, to leading order, the parameters of the quadratic part 



of the action do not flow, i.e. , 

Zn+l ~ Z n ~ Zq 

v Fn+ i « v Fn « v F0 (24) 

A n+ i S/i n+1 « 0(g) . 



In order to analyze the flow of the quartic coupling constants g^ a ' (ui) , we 
classify them in terms of qualitatively different channels. Besides the spin indices, 
a, a' , they depend on the 4 discrete momenta lj = {uj[, U4} . At the n th 
iteration step, each of these unit vectors can take iV„ ~ A^ 1 different values 
ojW i i = i i ...,iV n . However, the four momenta u are not all independent, but 
must satisfy the momentum conservation u\ + u 2 = ^3 + w 4 , as required by 
translation invariance. Studying the geometry of momentum conservation, we 
can subdivide the set of coupling constants into three different channels. There 
is a qualitative difference between two and more than two dimensions. 

In d = 2, cf. Fig. 6a, given the two incoming momenta Ui, U2, with 001+0)2 7^ 0, 
there are exactly 2 possibilities to choose the outgoing momenta, either uj 3 = a3 2 
and 0)4 = uj\ or u5 3 = uj\ and U4 = u 2 ■ In the exceptional case where 00% + u 2 = 0, 
one is free to choose arbitrarily one of the A n ~ N n discrete values for a; 4 ; d} 3 is 
then determined as its antipode, U3 = —Q^ (i.e. , in this case, not Q\ and u 2 are 
independent momenta, but to\ and co^ ) . 
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Fig. 6a 



Fig.6b 



In higher dimensions (cf. Fig. 6b for d=3 ), given the two incoming momenta 
uj±,uj 2 , with uji+u 2 7^ , one has 0(X d l ~ 2 ) choices for the outgoing momenta. One 
of the outgoing momenta, e.g. , u 3 can be chosen arbitrarily on the d — 2 dimen- 
sional sphere S d ~ 2 (u)i, u) 2 ) , £4 being determined as its antipode on S d ~ 2 {uj\ 1 uj 2 ). 
Again, in the exceptional case, where Q\ + u) 2 = , there are 0(A^ _1 ) choices, as 
in 2 dimensions ( i.e., the number of choices of outgoing momenta is larger than 
for a generic configuration by a factor of ~ A„) . 



Thus, in two dimensions, the couplings can be classified as follows : 
# d (o;i • u) 2 ,a ■ a') := g aa ' (ui, uj 2 , UJ3 = &2, ^4 = £1) 



g e (u x ■ u) 2 ,a ■ a') 



(25) 



9" 



(Wl , U 2 = -Ui , UJ 3 = -U4 , CO4} , 



cf. Fig. 7 . Here, "d" stands for "direct", "e" for "exchange" and "c" for "Cooper"; 
we use the convention that a = ±1 for "spin up" or "spin down", resp. . 



0J\ a 



a' ^2 



uj 2 a 



cr' Wl 



004 a 



a' -U4 



N 



UJ 2 



N 
Fig. 7 
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N 



U) 2 



21 



Note that because of rotational invariance the coupling constants, g(u), really 
only depend on the scalar product, u • uj' , of two vectors on the unit sphere, 
i.e., on the angle, u'), between them. Analogously, they only depend on the 
relative orientations, a ■ a' , (i.e. parallel or antiparallel) of the spin degrees of 
freedom. 

For an electron system in more than 2 dimensions, momentum conservation allows 
more independent coupling constants than the ones listed in eq . (^5|) . However, 
the analysis of their renormalization flow in leading order in j- is similar to the 
one for two-dimensional systems. To leading order in they neither flow nor 
do they influence the flow of the ones listed in eq . ([2~5|) . 

The contributions to the functional T(u, V n ) of leading order in can be 
determined by applying similar geometric considerations. However, now one has 
to study the geometry of momentum conservation for vectors lying in thin shells 
around the Fermi surface . 

In Fig. 8 , the contributions of leading order in y- are displayed (except that self- 
energy contributions to the inner propagator lines are omitted). A summation, 
XX') > over k> £ is associated with every electron loop ( for graphical rea- 
sons, the spin indices are omitted in the loops; we sum over all compatible spin 
configurations). 
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Fig. 8 



The renormalization of the coupling constants to leading order in g is deter- 
mined - for each channel - by the first two diagrams. Diagrams a) renormalize 
the direct and exchange channel for electrons with parallel spins and the direct 
channel for electrons with opposite spins, diagrams b) renormalize the exchange 
channel for electrons with opposite spins and diagrams c) the Cooper channel. 
The contributions corresponding to diagrams a) and b) only lead to a weak renor- 
malization of the direct and exchange channel and, solving eq.(^), one can show 
that - for sufficiently small initial conditions - the running coupling constants 
corresponding to these channels stay bounded. (This result remains true if one 
includes the leading order contributions of the higher vertices shown in Fig. 3 , 
cf pi) . 

The renormalization of the Cooper channel is - to leading order in -p- and omit- 
ting self energy renormalizations of inner propagator lines - independent of the 
direct and exchange channel : 



An 



(26) 



l + (3gC(£;a-a') 

where (3 is a positive constant of order unity, and I = 0,1,2,... refers to the 
angular decomposition of the Cooper channel : e.g., in d — 2, the trigonometric 
decomposition 

oo 

gS(u 1 -u 2 = cos<&) = Y.&)cos{B) 

e=o 

is used ( in d = 3, the Cooper channel coupling function is expanded in a sum of 
Legendre functions, rather than trigonometric functions ). 
The solution of eq. (|26"D is 

~ i + in + WtfV;*-*) • (27) 

Suppose that one of the Cooper channel couplings, e.g. Qq^*) , is negative, 
(£* = 0, 1,2, ...) . Then we encounter a singularity after n* = a \ g o^»^ iteration 
steps, corresponding to an energy scale ^F^{j^) n * ■ Thus, perturbation theory 
breaks down when n — > n* . This is the celebrated superconducting instability 
(for a more detailed discussion, see |17|). 



However, one should note that this instability only appears if the Fermi surface is 
symmetric under reflections at the origin, i.e., keS^' 1 implies that -k € S^ 1 . 
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If the reflected Fermi surface intersects the original Fermi surface transversally, 
the superconducting channel is not renormalized to leading order in y-. 
If all initial values of the superconducting couplings are positive, i.e., g$ {€) > 
, £ = 0, 1, 2, ... , then, following eq. (|27|) , they all appear to flow to zero. How- 
ever, in general, higher order corrections (in g and y-) to eq.(|26|) invalidate this 
conclusion; this is the " Kohn-Luttinger" - effect |18| . An analysis of this effect 
can be found in fl|, |T4| . 



In the following, we assume that the Cooper channel is turned off, as e.g. in 
systems of electrons coupled to magnetic impurities with broken parity- and time 
reversal invariance fl[] . Then, for short-range interactions, all remaining coupling 
constants are marginal, i.e., g n +i ~ 9n- (This result turns out to be stable against 
adding higher order corrections). In this the coupling constants in the 

effective action S^ ff are suppressed by a factor , our calculations indicate that 
the system flows to a Landau Liquid Fixed Point, as n — > oo. In the following 
sections we confirm this expectation by calculating the electron propagator using 
bosonization. The effective action S^ ff obtained by the RG analysis serves as an 
input for this calculation. 

For long-range interactions, the situation is more complicated. For long-range 
density-density interactions of the form displayed in eq.(^), the effective action 
Sq S at a large, initial scale A has the form given in eqs. (|I2"D - flT5| ) . The ef- 
fective interaction potential can be replaced by a set of coupling functions, as 
described in eq.([T7|). Due to the singularity of the interaction potential V"(|p|) 
at \p\ = 0, the interaction processes with u} 4 ^ uj\ are less important than the 
direct scattering processes, where = Q\ and uj^ = Q2 ■ In this channel, the 
singularity contributes a supplemantary factor Aq to the quartic term ( |T5D in the 
action which can partially or completely compensate the factor —x=t ■ However, 
resumming all diagrams of leading order in leads to " screening" which renders 
the long-range interaction effectively short-ranged. Therefore similar results as 
for the short-range case are expected to hold. We shall confirm this expectation 
in the approximation obtained by bosonizing the system. 

For long-range transverse current-current interactions, as they occur in Quantum 
Hall fluids at filling factors v = |, |, ... , the screening mechanism is ineffective. 
Calculating the electron propagator by the bosonization technique, we shall ob- 
serve the possibility for a deviation from Landau Liquid behaviour, depending on 
the exponent a which characterizes the singularity of the interaction potential in 
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momentum space : for a > d — 1, we argue that the system is a MFL (similar 
predictions have previously been made in 0, f|, [5], ||). 



3 Effective Gauge Field Action and Bosoniza- 
tion 



It is easiest to understand the meaning and accuracy of "bosonization" by calcu- 
lating the scaling limit of the effective gauge field action, W(A) , where A is an 
external electromagnetic vector potential. From the effective gauge field action 
one can determine the (connected) Green functions of currents by differentiat- 
ing with respect to the gauge field A . Calculating the Green functions for the 
electron fields , ^ is more complicated and is accomplished in the next section. 

The external gauge field A p , p — 0, 1, d , is coupled to the electron system 
by replacing derivatives in the free action S° (\I/*, \1/ ; /x) , eq.(|TJ), by covariant ones 
D p (A) := dp — ieAp ("minimal coupling"). Here e is the elementary electric 
charge (we choose units such that % — c — 1) . Eq. ([[]) is then replaced by 

S° A) = S (^*,^;fi) + S J {V*,V;A) , (28) 

where 

S J (^,V;A) = fd d+1 xJ2A p (x)f(^,^,A)(x) , 

and the current density, j p , is defined (in euclidean space-time) by 
j° (x) = -ie^*(x)^(x) 

(29) 

A) (*) = ±[*'(x)D t (A)*(x)-(*(x)D t (A)y*(x)) , 
for £ = 1, d . 

The effective gauge field action W(A) is obtained by integrating out the degrees 
of freedom of the electrons, 

W(A) := -log j^ 1 J e -[s°(*%* ; *A) + ^(*-,*)]j ; (30) 

with S (\I/*, \1/) given by eq.(Q) . It is the generating functional for the connected 
Green functions. At non-coinciding arguments, one has that 

n r In \ con 

n 7i ? v w(a) = (-lrmriv^Aixi)) . (31) 

i=1 0Ji pi {Xi) A=Q \ i=1 / A=0 
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First, we consider a system of non- interacting electrons. The scaling limit, 
W°(A) , of the effective gauge field action W°(A) has been calculated in |L5 . 
Here, we just sketch the essential ideas and recall the main results. 
We expand the effective action W°(^4) - eq.(^), with 5=0- in powers of 
the field A : 

oo -i „ n 

W°{A) = E - / U dx t +1 ('"' " (•<•. '•><) A Pl ( Xl ) ■ ■ ■ A Pn (x n ) . (32) 

n=l ^" i=l 

The expansion coefficients, C , are given - at non-coinciding arguments - by the 



current Green functions, cf. eq. fl3~T|) . 

Next, we map the physical system in a space-time region to a reference system 
in the region A^ = A , where A > 1 is a scale parameter, and /i is kept fixed. We 
shall be interested in the asymptotics when A — > oo. Under the rescaling map, 
points in A^ transform as 



x 



A£ e A {x) i— > £ e A 



The gauge field A. x '(x) - which probes the response of the electron system to 
small external electromagnetic fields - is chosen as follows : 

Af\x) = i a ,(0 , (33) 

where a p (^) is an arbitrary, but fixed function on A. Thus, the gauge field scales 
like the momentum operator. In the following calculations we consider the for- 
mal thermodynamic limit A — > R d+1 . To construct the scaling limit W°(a) of 
W°(A^) , we study the asymptotic form of the current Green functions, eq.(pl|), 
in the limit A — > oo , using eq. (§) , and plug the result into eq. ([52]) . We define the 
scaling limit VVf(a) as the coefficient of the most divergent term in an expansion 
of W°(A W ) in powers of A and A -1 . 

In the scaling limit, one encounters ultraviolet divergencies in the perturbation 
series of W°(a) . To fix a resulting ambiguity, we use standard Ward identities 
implied by gauge invariance, i.e., 

d 

r<Pl,-,Pi,-,Pn (- ™. „ \ _ n 

w -"n •••) ^n) w 



dx 



The result of our analysis can be summarized as follows. In a calculation of 
electron Green functions at distance- and time scales of order A, the free fermion 
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action S° fi?*, fi) can be replaced by the following approximate action : 

=: S°({A>) ■ (34) 

The symbol indicates that the approximate action reproduces the large 
distance- and time asymptotics of electron Green functions to leading order in j. 
The sum, , extends over the discrete set of unit vectors, uJj, i — 1, N ~ 

A^ 1 , whose endpoints lie on the surface, <Sf -1 , of the d— dimensional unit sphere 
(cf. Section 2). The Fourier modes V>£(fc) have support on I u := R x B w . In 
the second line of eq.(|34|), we introduce an (arbitrary) partitioning of iSf -1 into a 
positive, S + , and a negative, 5~ , hemisphere and denote the ray through u , for 
uj G S + , by [a;] := {cD, —a;} . This allows us to use the "relativistic" notation 

*M : = ( ^ ) *M : = *M7° = ( C ) , (35) 

with 7° = <j\ and 7 1 = o<i . In the following , we assume that electron spins are 
frozen in a fixed direction (we shall omit spin indices) . 

In determining the effective gauge field action W°, the term S J (^*, ^ A) , 
coupling the gauge field Ay\x) = j a p (£) to the electron fields, can be replaced 
by 

S J A^) « E H«o(()-^)]' 

• E eW-*'* : : 

=: S J (tyJ};a) . (36) 

Again, the symbol "~" indicates that, in a calculation of electron Green functions 
at distance- and time scales of order A, this approximation yields the leading 
contribution in j . As mentioned above, in the expansion of the effective gauge 
field action W° in powers of A, we evaluate the expansion coefficients - i.e. the 
current Green functions - only at non-coinciding points. The omission of self 
contractions is indicated by the normal ordering of the product of electron fields in 
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eq.(pq). The local terms of the expansion coefficients are determined by requiring 
gauge invar iance. 

Given a ray [u] , we define a current density 

JHiO ■= E {e ikF ^-^ : VMOlMO : + e ik ^~^ : ^(O^-AO ■} , 

(37) 

corresponding to electron motion in the direction of —u , and a current density 

7[«](0 : = E {e-*'^^ + e"*'^^ : :} 

(38) 

corresponding to electron motion in the direction of u . 

Using eq.flrof), one can establish that the quasi 1 + 1 dimensional current densities 



jUO ■= I (iM(0+j M (0) 



(39) 



for [to] E [Si x ] , are related to the d + 1 dimensional current density, j p , p 
0, ...,d , defined in eq . (p9|) , by 



j°(x) := j°(**,¥)(s) « A- d E jUO 

(40) 

j\x) := A) (x) « X- d E , 

A=0 ues+ 



where £ = 1, ...,d . Eq. ( fjO|) holds for current Green functions to leading order 
in i at non-coinciding arguments . For simplicity, the electron charge, — e , has 
been set to —1. 

Eq. (|36|) can be written in the suggestive form 

S J (**,*;AW) « 5 J ({A};a) := E E / ^ , (41) 

£=0,1 

with 

a o (0 := a o(0 , <(0 := v F ua(£) . 



One observes that - to leading order in the inverse scale parameter A -1 - the 
original d + 1 dimensional system decomposes into independent, quasi 1 + 1 di- 
mensional subsystems, one along each direction [u] of R d . The 1 + 1 dimensional 
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subsystems describe relativistic fermions moving along the direction [u>] with ve- 
locity rkvp. Hence, the calculation ofW®(a) is reduced to the calculation of the 
effective gauge field action of a family of independent Schwinger models . 
One obtains 



where we used 



1 1 ,k 



Fxd- 



\ d 1 7T 2ir vf 



E / dd+1 t E «2 T (O«2 T (0 



£=0,1 



2 JnxB ~_ n 



p,a=0 



(42) 



(0 == E [Sb, 



c 



U B°C 



C=0,1 



«c(0 



and 



(43) 



Note the non-trivial (and important) fact that W° is only quadratic in the gauge 
field. 

Invariance under gauge transformations, space rotations and parity - or time - 
reversal implies that the (euclidean) polarization tensor Ilf 7 has the general form 



k, 



Uf(k) = Uf(k) = -^ui(k) 



(44) 



n?(fc) = nl(k)U j -^)+nt(k)^f i,j = i,..,d 



The calculation of the two independent functions Ilf and II*, for non- interacting 
electrons, yields the result : 



ni(k) 



nt(k) 



^k F v F 

7T 1 Vp fc2 



'-k F v F 



(1- 



fc 



-Hi? A; 



vpk 
k 



\ko\ 



yJkl+(v F k) 2 



arctg 

Ui(k) = e -k F v F 



d = 2 
d = 3 



(45) 



arctg 



-n'(fc 



d = 2 
d = 3 
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Comparing this result to the small-fc asymptotics of formulas found in standard 
textbooks, one finds agreement - up to the " diamagnetic term" of the transversal 
part which is proportional to \k\ 2 and yields lower order corrections in j. 

Eq.([E|) is the result of resumming the leading contributions in the (formal) 
expansion of the effective gauge field action W°(v4) in powers of k- 
By power counting, the possible diagrams contributing to leading order in j to 
the effective gauge field action are the following ones : 







+ (3 + 



Fig. 9 



Straight lines stand for the electron propagators 

1 



Gl{k) 



ik — vpujk 



and crosses for the gauge field insertions (note that, in contrast to the diagram- 
mar introduced in the context of the RG, the electron propagators are understood 
to be integrated over the entire domain I u = R x B w ) . 

However, eq.(fl^) states that, because of the special form of the propagators G°, 
only the first term in Fig. 9 with two gauge field insertions contribute to W°(a) . 
In addition, gauge invariance requires the introduction of local Schwinger terms 
which couple, along a ray [u] , right- and left-movers. 



Remembering that, in the 1 + 1 dimensional Schwinger model, the calculation 
of current Green functions in the fermionic theory can be reproduced by express- 
ing the current densities in terms of a free, massless bose field (see, e.g., |L5j, [19 
and references therein), we introduce following identifications : 



2 k F ti 

7^ 2 d ^'(0 
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(46) 



.7M (0 



2 ,k F d-i-^jj 



0r 27r' 

for a; e 5 + , with [a;] = {to, —uj} denoting the corresponding ray, and 



^ 1 f . d ^ 



v 



I d _ d 

2{'^ + VF0J Ti 



Eq.(|4"6"D is equivalent to 



% , ki 



2 e 



BC 



y^i" 27r' 
for 5, C = 0, 1 and u e S + . 

The action of the set {</?m(£)} °f bose fields is given by 



S°(W H }) = ~E / - <? d+1 ^M(-fc) k 2 + (^) 2 | M (*) , (47) 



where the Fourier modes, (p^k), of the bose field <^m(0 fi ave support on Rx B^. 
( For a given value fco, the modes 0[ u ](ko, ■) have a compact support. We choose 
this support to be B w , but, in principle, there is no need to choose a support 
identical to the one of the fermion modes ^(ko, •). ) 



One can verify that 

s;'/p{ w }e- s '«"H}) jn 



(48 



in the sense of distributions. 



The representation (£5^) of the one-dimensional current densities, , in terms 
of the bose fields, ip^] , is implied by the local conservation of the electron cur- 
rent, eq. (|40|) , along each ray [u] seperately. These formulas are equivalent to 
Luther-Haldane bosonization [JlTJ, |l^, 

The identification ( j4^) reproduces only the leading term of the fermionic pertur- 
bation theory, i.e. : 

W°(a) = " l^''-%{^ 1 /p{^ M }e-[ s0 ^N})+^toMha)]} , (49) 
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with 

/ <r r 'X o^i):'-" — ( /w 



The symbol " lim^oo " stands for determining the coefficient of the leading term 
of the expansion in A and A -1 . 

The derivatives d^(p\ w ] and d ip^] of the bosonic fields ipy] describe the current 
density fluctuations, j^j and jr^i , determining the scaling limit of the system. For 
a fixed direction [u] they are composed of electron modes with momenta near 
the points kpQ and —kptd, resp., on the Fermi surface. In the scaling limit, they 
probe the Fermi surface only locally around the points kpoo and —kpu . 

This implementation of bosonization is a special realization of a more general 
formalism, presented in ]20| . 

Abelian gauge invariance implies the local conservation law 

d p f(x) = . (50) 

Instead of expressing the current density in terms of the elementary fields \&* , $ 
and imposing the constraint ([50]), one can introduce new field variables which 
guarantee eq. (|50|) by construction. For d + 1 dimensional currents, j p , one needs, 
in general, antisymmetric gauge forms of rank d — 1 . However, considering only 
the scaling limit of the system, we can take advantage of a substantial simplifica- 
tion : the fermionic theory decomposes into a family of quasi 1 + 1 dimensional 
subsystems ("dimensional reduction"), one along each direction [a;] of R . Each 
subsystem describes quasi 1 + 1 dimensional, " relativistic" electrons. Gauge inve- 
riance has to be fullfilled for each subsystem - i.e. in each direction - seperately, 
and implies the conservation laws for the associated currents, jy , B = 0, 1 , 
i.e. 

+ = . (51) 

Expressing the currents jy, in terms of the free, massless bose fields (py guaran- 
tees that eq.(|5"TD holds. 

Next, we study the effective gauge field action for interacting systems. 
In section 2, by using the RG group method, we successively eliminated the elec- 
tron modes with momenta outside the shell £l\ n :={pe R d ; p — fcifjC < jf- | 
of thickness ^ fop around the Fermi surface, in order to determine the effec- 
tive action S n for the remaining modes with momenta inside the shell fl\ n . The 
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resulting effective action S n describes the physics at energy scales smaller than 
vf^~- Under the assumptions specified in Sect. 2, one can reach very small energy 
scales iff -p <C VpkF before the form of the corresponding effective action S n 
differs considerably from the original action S. 

In this section, we consider "spinless" electrons and suppose that, in addition, the 
Cooper channel is turned off. Then the terms S% and S* , quadratic and quartic 
in the fields ^(k), are given by 

~ E / ^CW^ (ik Q -v Fn uk) Mk) (52) 



and 



9 \d-l 72 T 
^ /v n , , , , c-od-1 



with 



feC^^C*®)^^)^^) , (53) 



# n (^ 2 ;i) = {^(wiwa;^) - c^(^ 2 ; £)} . (54) 



The "direct"- and "exchange" -coupling functions g d and g e have been defined 
in eq . (p5|) . The symbol "~" stands for "equal to leading order in j-" . 
Contributions to S n corresponding to higher orders in the Taylor expansion of 
the coefficient functions of the quadratic and quartic terms in the momentum 
variables, or contributions involving more than 4 electron fields are neglected. 
By engineering scaling they are irrelevant, cf. Sect. 2 . For systems with a short 
range two-body interaction potential V, the coupling functions g n (u}iuj2', x~) can 
be replaced - to leading order in y — by the coupling constants g n (u;iu; 2 ) ■ We 
first restrict our attention to such systems. 

The term S^^ip^}; a) coupling the electron fields to the gauge field a has the 
form 



S^lha) « £ J2 K +1 £a£(£)j&(£) , (55) 

o, ei S+ 5=0,1 

where (aft, a") := (ao,VF n 3 • a) , and the quasi 1 + 1 dimensional currents 
are defined as in the non- interacting system, cf. eqs. (|37| ) - (0), but with 
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the products VC'(£)V^(0 of the electron fields replaced by ^ ^jy (0^w(0 ( m 
the non-interacting system Z n = 1). A Ward identity relates the renormalization 
of the vertex ([55]) to the renormalization of the electron propagator, preventing 
the appearence of a new, independent renormalization factor. 
The quartic interaction term given by eq . (|53|) can be expressed in terms of the 
quasi 1 + 1 dimensional currents : 



g n [-LOiU 2 ) 



Jlvi] (£) Jfa] (V) + 3[ U1 ] (0 j[wa] fa) 



(56) 



with 



JR ^|fe|<-f- 



Replacing the quadratic term of the action by the bosonic action S°({ip[ u ]}) 
given in eq.([4~7D and inserting the bosonization identities (|46|) for the current 
densities in eqs. ( |55| ) and (|56|), we end up with a gaussian bosonic theory The 
bosonized versions of eqs. (|55|) and ( |56"D are given - in momentum space - by 

Si({<PMh") = E E / a d+1 ka%(-k)e BC ^(^)^k» om (k) (57) 

W&S + B,C=0,1 V 71 " 

and 

1 , 1 



X ({<?»}) = w & EL«' 



with 



+ g n {uSii2 2 ) n^](~ k ) (k^k^ +r i r 2 ) V[W2] {k) 

(5J 

1 - — w 1 

/c w = -(— ifc + v F ujk) and = -(iA; + VpCuk) 



With the aim of calculating the scaling limit, Wj / (a), of the effective gauge 
field action of the interacting electron system, we replace the fermionic action 
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S n , given by eqs. (p2[) - (0), by their bosonic version, given by eqs. (|47|), (|57D 
and (p8]), calculate the corresponding effective gauge field action, W y (a; V n ), and 
determine the coefficient, Wj / (a; "P n ), of the leading contribution in an expansion 
of W v (a;V n ) in powers of A n and A" 1 , as \ n — > oo. V n stands for the set of 
parameters characterizing the fermionic action S n . As the bosonic theory is gaus- 
sian, it is easier to study than the fermionic theory where the interaction term is 
quartic in the fields ify. 

But does the action Wj / (a; P n ) derived by the bosonization procedure repro- 
duce the correct result that one would obtain by iterating the RG transformations 
and by determining the resulting fixed point action (n — > oo) ? 
In order to answer this question, we have to clarify the approximations involved 
by bosonizing the electron system described by the action S n . The approxima- 
tions can be characterized in terms of the formal perturbation expansion of the 
fermionic theory in powers of A" 1 . In fact, by bosonizing an electron system 
with action S n , we resum all leading order contributions of the fermionic per- 
turbation expansion of the effective gauge field action in powers of A~ x - except 
self-energy renormalizations of inner propagator lines. The diagrams reproduced 
after bosonization are the following ones : 




Fig. 10 

As above, straight lines stand for electron propagators and crosses for gauge 
field insertions. The interaction -frr g n (uiU2) vertices are represented by dots. 
In the diagrams in Fig. 10 , each factor per interaction vertex is compensated 
by a bubble-summation, XX') • 

These diagrams correspond to the so called RPA approximation . 

In principle, to leading order in a (formal) expansion in j-, the propagator 
lines in the polarization bubbles in Fig. 10 are renormalized by "cactus diagrams" 
of the form : 
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Fig. 11 



For systems with a short-range two-body potential V, their contribution is - to 
leading order in y — a constant, i.e., they only lead to a displacement of the 

-An 

chemical potential. As the polarization bubbles only depend on the difference of 
the momenta of the two inner propagator lines, they are not modified by such a 
displacement of the chemical potential. 

Does this result imply that, by taking functional derivatives of W% \a;V n ) 
with respect to the gauge field a and by evaluating the resulting current Green 
functions, cf. eq.(|3l]), for large arguments — = A n |£^ — C \ > with 



X n — > oo, one obtains the leading contribution in j- of the current Green functions 
of the interacting system ? One must remember that we use the fermionic action 
S n as an input of bosonization. Because of the linearization of the pieces of the 
Fermi surface contained in the boxes B w , this action describes the properties of 
the system only correctly at momentum scales between ^ and fe, as argued in 

An A n 

Sect. 2 . 

It follows that the current Green functions derived from Wj / (a; V n ) reproduce the 
leading order contribution in y - of the current Green functions of the interacting 
system for arguments \x^ 1 ' —x^>\ = A n |£^ — C | of order between A n and A^. In 
order to explore the current Green functions of the interacting system at larger 
distance- and time scales, we first have to to iterate the RG transformations 
further, with the aim of deriving effective actions S m , m > n, describing the 
properties of the system at larger distance- and time scales. Such an effective 
action S m can serve as a new input for the bosonization procedure. 

In Sect. 2, we have analyzed the flow of the set V n of parameters characterizing 
the effective action S n for electron systems with a short-range two-body potential 
V. All the parameters of the set V n - except the scale factor A„ - tend to finite 
values V* in the limit A n — > oo. This implies that 
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W;(a) = " lim " W v (a;P n ) , (59) 

n— >oo 

i.e. the scaling limit of the effective gauge field action is given by 
W?(a) = " lim" -log I-' 1 [ V{<p { u] } e -(^M})+S* J «*>]};»)) e -^({^,}) 



(60) 

where the action is defined by eqs. (|47l), fl57j) and (|58D . Again, "liniA n ^oo" stands 
for determining the coefficient of the leading term of the Laurent expansion in A n 
and A^ 1 . 

Carrying out the Gaussian functional integral, we find for the leading term pro- 
portional to A^ -1 : 

w» = \ fd d+1 k £ a p (-fc)n&(fc)a.(fc) • (61) 

p,CT=0 

The polarization tensor Hy* (k) has the general form described in eq . (fill) , where 
Il V:ls (k) is given by n*(&), and Tly^(k) can be expressed as a Neumann series 
in g*(uu'). The explicit expression for the Neumann series is rather complicated 
in the general case (but straightforward to determine). We only write down the 
explicit expression for YI Vjt (k) in the (somewhat artificial) case when ^(^1^2) = 
g*. Then one obtains 



1 + k i^mk) 

(The functions U 1 / are defined in eq.(fl6])) . 

For systems with a long-range two-body potential V, the renormalization flow 
of the effective action S n is not yet fully understood. Given the effective action S n , 
we can calculate calculate the corresponding effective gauge field action W(a, V n ) 
by bosonizing the system. This yields a result for the current Green functions 
in the domain of validity discussed above. In order to derive results at larger 
distance- and time-scales, one is obliged to make assumptions on the flow of the 
effective action S n , as n — > 00. 

In the following, we discuss the calculation of W(a, V n ) for three examples of 
systems with singular interactions : 
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i) Long-range, density-density interactions 



We consider interactions of the type shown in eq.(f2|) with a long-range inter- 
action potential gV, i.e., one whose Fourier transform, 0VY[p|) = -Jhr ^S- , 

k F P 

becomes singular at \p\ — . The exponent a and the coupling constant g are 
supposed to be positive. We assume that the effective action S n at an energy 
scale vpn^- has the form specified in eqs. (|52|) - (|55|) and that all the parameters 



of the set V n characterizing the action S n - except the scale parameter X n - tend 
to finite values in the limit n — > oo. In particular, the set of coupling functions 
flwG^j t~£) * s supposed to be related to the initial interaction potential by eq. (|I~7|) , 
for arbitrary n. We can neglect the exchange channel with respect to the direct 
channel. In the direct channel - contrary to the case of short-range interactions 
- we have to retain the (singular) dependence on the small momenta, -j^k , as 
indicated in eq. (|T7D . Hence, in eq. Qolf) , the dominant coupling constants are given 
by 

k F \ n 



g n (^2]^k) = k d F 1 g n V{^\k\) = g n 
By bosonizing the system, one obtains the result 



k 



(62) 



W D {a-V n ) = ^X^Jd^k J2 a p (-k)U^(k;V n )a a {k) , (63) 



p,CT=0 



where the polarization tensor YI p ^(k;V n ) has the same form as in eq. (|44|) , with 
the two independent functions Yl l D and LT^, given by 



IlUk-V n ) = W(k) =- , — = 

1 + %9nV{\±\)W n {k) 



(The functions Ilff are defined in eq.([46|), with vp — > vfu-) 

As expected, this result coincides with the result of an RPA calculation 



ii) Tomographic Luttinger Liquid 

We introduce interactions which describe singular direct (or "forward") scat- 
tering processes, where the quasi-electron with momentum near kpti interacts 
only with one with momentum near —kpuJ ■ This leads to the coupling constants 

n (tfitf 2 ) = Xt 1 (Svr)^ 1 g n . (64) 
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By inserting eq. (|6lD into eq. (|BD[) , one obtains 



p,cr=0 

The polarization tensor 11^ is renormalized by an overall factor (l + f^) with 
respect to the non-interacting system : 

n&(*;P„) = (i + f£) _1 iC(*) ■ (65) 



iii) Long-range, transverse, current-current interactions 

In the original system, interactions between transverse currents lead to an 
additional term in the action of the form 

S T (^,V) = -I / d d+1 x I d d+1 y^j^x) 9V{\x-y\)5{x Q -y Q )j k T {y) , (66) 
1 J J k=i 

with 

1=1 \ l^i=0 °i J 

A physical realization of such a system is a QH system at filling factors v = 
2,|,..., which can be described as a system of free (composite) fermions interact- 
ing via long-range, transverse, current-current interactions as diplayed in eq.(|66|). 

cf. ©si . 

One determines the effective interaction for the low-energy modes whose momenta 
lie in a thin shell Q n around the Fermi surface following the same procedure as 
for longitudinal interactions. Using eq. (f41|) , one obtains 

sum) = -f^r e / T a d+i kv(\^-\)- 

■ E ^\3U(-k)Vl{k)uijl 2] {k) , 

i,j=l 

(67) 
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with 



k{kj 

3 k 2 



and 



V{\P\) 



t.d-1 
F 



P 



For a QH system with unscreened Coulomb interactions, the exponent a is equal 
to 1 ; if one assumes the Coulomb interactions to be screened, the exponent a 
turns out to be 2. 

We propose to calculate the effective gauge field action of this system by replacing 
the fermionic currents, ^(V^j VVl) > ^ their bosonic versions, J^](<£[tj]) , given 
in eq.(^6|). We have to make assumptions analoguous to the ones stated in (i) . 
One obtains 

d 

d d+1 k Y, CL p {-k)Tl^{k-V n )a a {k) , (68) 



W 1 (a;V n ) 



1 \d-l 

2 K 



p,(T=0 



with 



II e T {h,Vn 
n^(/s; V n 



K{k) 



K(k) 



1 



1 + W n (k)g n V(\l\) 
Again, this reproduces the result of an RPA calculation. 



4 Electron Propagator 

In this section, we determine the bosonic expressions for the electron fields 
in order to calculate the electron propagator for interacting systems. We shall 
bosonize each one of the iV ~ A d_1 component-fields ip^ seperately (/0" stands 
for ip* or if) ). More precisely, for each ray [u] = {Co, — u} , we express the 
pair ^Ej^-w i n terms of a bose field (p^. This is accomplished by applying the 
well-known bosonization formalism for 1 + 1 dimensional relativistic fermions 
summarized in appendix A . However, one has to cope with a subtlety arising 
from the dependence of the quasi 1 + 1 dimensional electron fields ipt on the 
components of the momentum perpendicular to the direction [u] . 

We start our discussion with the non-interacting system for which 

-(&«(*) fco(A0)° = i^) d+1 5 {d+1) (k-k>) ljj (k) . 

tk — vpuj (l 'k ^ ' 

(69) 
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In comparison to eg. (|Tl"D , we replaced the boxes J5 W (») by cubes, Q w (i) , with sides 
of length kp ( the exact shape of the integration domain is irrelevant) . Because 
of the isotropy of the electron system, one observes a non-trivial dependence on 
the components of k only in the 0— and a;— direction. Propagation takes place in 
the radial direction, u. This suggests a decomposition of the quasi particle fields 
■^£(0 i n t° tensor products 

vico =xM(a)®^i[|(Mi) , (70) 

for G 5f- x . 

The "radial" fields, ip^ , describe 1 + 1 dimensional, relativistic electrons, whereas 
the bosonic "angular" fields, xm(C-J-) > J us t guarantee momentum conservation in 
the perpendicular direction. 

The angular fields Xm(£l) are Gaussian, of mean 0, and their propagator (co- 
variance) is given , in momentum space, by 

(x M (fc±)x K ](^)) ± = S m (27r) d - 1 5^(k ± -k' ± ) • • (71) 

Their Fourier modes, Xm(£l) > have support in \ 

The radial fields can be bosonized by applying the standard 1 + 1 dimensional 
formalism . For each ray [u], one uses the following identifications : 

V>u4£o,£|[) : = (Co, < — > 
€i(£o,£||) := C||(6,en) — 
^2(6, := ^||(6o,£||) < > 

c 2 (e ,en) := Cii(6,en) «— > 

where D[ u ] is a disorder operator and the normal ordered exponential of is a 
"vertex operator" . The precise definitions of the expressions on the r.s. of eq . (ff2|) 
appear in appendix A . 

The fields ^L](£o>£||) are 1 + 1 dimensional free, massless Bose fields with an 
action given by 

^iML]) = \ / Rx[ _^ ^ *Mfcii ^h(-^o, -*«) 0o + (^i) 2 ) ^L#o, . 

(73) 



£>m(£o,£ii; 1) : 
1 D [w] (£ ,£|i;-i) 

v^](eo,en;i) : 

K) 4 

1 D [w] (e ,eii;-i) 

7T H 



e *v^V'L](£o,?n) . 

g-iv^JLj^O'Sll) . 

. e *>/^[L]«0.«l|) ; 
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The Fourier modes ^[L(fco, k\\) have support on R x [— , ] . 
One can verify that from the bosonization formulas, eqs. (|7(])-([72|), for the elec- 
tron fields the ones for the current densities, eg. (flop , follow. In appendix B, we 
show this in an example. 

It is important to understand the relation of the radial bose fields wL(£oj £||) 
to the fields <P[ W ](£) introduced in eqs.(|46|)-([47l), which are related to the fermionic 
currents. 

Following eq. (|47D , the action S (<p[u]) of the bose field tf[ u ] is given, in momentum 
space, by 

S°(<PM) = ~\ \ - a d+1 k^ H {-k) (kl + {v F ktf) <p[u](k) , (74) 

i.e., the field ^r^i propagates only along the direction [u] . It follows that the 
propagators of the ip^ -fields are related to the ones of the </?L- fields by 

(v>](6,6||,<?-l) n^'](Vo,V\\,V±)) = 

s W] c^cl - v±) ( ^[L](Co, e«) vUivo, m) )" . (75) 

with 




d d-i k e ik ± u 



Remark : Comparing eq.(|T4j) to eq.(|7^), one can decompose the field 

into a tensor product y>[ w ](f) = <£[L](£o, £||) ® , where Xm(H) 

coincides with the bosonic Gaussian field introduced in eqs. (ffOl) and 

(0)- 

Next, we study the effect of interactions of the form displayed in eq.(|SB|) . 
As discussed in the preceeding sections, we use the effective actions S n given in 
eqs. Q52] ) and as an input for the bosonization calculations. The action S n is 
the result of an RG analysis where the electron modes with momenta p outside 
the shell fl n have been integrated out. It describes the properties of the electron 
system at energy scales smaller than VFn^r- Under the assumptions specified in 
Sect. 2, one can reach very small energy scales VFn^- -C VFnkF before the form of 
S n deviates from the one given in eqs. (|52| ) and (0). The (remaining) fermionic 
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degrees of freedom are described by the fields W. These fields ^"(A^) can be 
decomposed into N n ~ X d ~ l independent components : 

(76) 

where the Fourier modes ifiu(k) have support in Q w (cf. eq.(|H])) . Momentum 
conservation guarantees that (ipuiip^') v ~ <L,o/ , so that the interacting propaga- 
tor, G , splits into the contributions of the component fields tfj^ : 

G(A n (£ - t?)) » E e*^^"* A; d G w ,„(£ - ry) , (77) 

where G u>n {^ — rj) are the propagators of the component fields calculated by 
using the effective action S n . The symbol indicates that the equation holds 
to leading order in an expansion in y-. 

As in eq.(|[), this formula should be regarded as a discrete approximation to a 
continuous angular decomposition of the interacting propagator, 

G(A n (£ - ry)) « / d*" 1 * e *^(f-« A; 1 Gj> fB - ifeMf- v)) , 

for large arguments A n (£— rj) , where the radial propagators G^ n (^o — Vo, — V) 
only depend on the 0-component and the component parallel to the direction uj. 
Hence, for the interacting system, the decomposition ( |7D| ) of each component ipf^ 
into a tensor product of a radial field (£o, £||) and an angular field xm(£i) - 
with the propagator defined in eq.flTTD - remains valid. 

We calculate the propagator G^ n (£o, £y) in the radial direction by bosonizing the 
system, i.e., we use the identifications Q72|) for the component fields and replace 
the fermionic action S n ({ipl,}) by its bosonic version <Sn({<£V]}) given in eqs. fl4"7| ) 
and Q58|). Then, the calculation of the propagator Crjj, n (£o>£||) is reduced to eval- 
uating an expactation value in the interacting bosonic ground state of a product 
of a disorder- and a vertex-operator in (plj, (see eq.([72])). In order to calculate 
this expectation value, we have to determine the action, S'n^L); °f the radial 
bose field from the action ^({^[w]}) of the family {ip[ u \ , [a?] G <S + } of d + 1 
dimensional bose fields. 

Through the interaction (0), the field ip^, for a given ray [uj], is coupled to all 
the other fields <£>[ w /], for [a;'] G 5 1 . Integrating out the fields <£>[ w '], for [a;'] ^ [u], 
we obtain the ("effective") action S n ((p[ u ]) of the field (p^ which is still Gaussian : 

VhW , (78) 



a d+1 k 



M (-k) Cfrfjz) 
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where 

= (0M(-*O0M(fc))n 

is the propagator of the Fourier mode <P[ w ]{k) in the interacting system. We 
determine the action S'nCvLi) °f the ra dial field ip 1 }-, (£ , by averaging the inverse 
propagator C'^(k) in eq.(ff^) over the components k± of k perpendicular to 
the direction [w] : 

1 [ dk o L kF _ kF &h <f>f u ](-fo, ~h) h U( k °> k \\) £[l]( fc o,fc||) , 

JR. J\ — ^ — , — o — ] 

(79) 



2 ' 2 



with 



-i 



Our procedure to calculate the radial, fermionic propagator GjJ, n is summarized 
in the following formula : 

Gfto.nCCo-^fll -»7[|) = -(Vw(Co,£[|)VCa(»to>»7||)}!l « 



(80) 

where a = 1, for —a;, and a = 2, for +a;. The disorder field D[ u ] in the bosonized 
expression of the electron field guarantees the correct anticommutation relations, 
regardless of the nature of interactions. 

Before diving into explicit calculations, we have to clarify the meaning of the 
symbol "~" in eq.flSOp, i.e., we have to clarify the approximations involved in 
calculating the electron propagator by means of bosonization. 
As in Sect. 3, where we have analysed the calculation of the effective gauge field 
action, we characterize our approximations in terms of the formal perturbation 
expansion of the fermionic theory in powers of A" 1 . By bosonizing the electron 
system with action S n , we resum all leading order contributions of the (formal) 
fermionic perturbation expansion of the electron propagator in powers of A,^ 1 - 
except for self-energy renormalizations of inner propagator lines (as in Sect. 3) . 
The diagrams contributing to the propagator G Wjn (fc) that are reproduced by 
bosonization can be found in the following way (we restrict our discussion to 1PI 
self-energy contributions) : 

First, draw all diagrams renormalizing the bare w-propagator line that are com- 
posed only of cu-propagators and contain loops of at most two propagator lines, 
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i.e. "bubbles". Then add all diagrams that are generated by replacing an oj- 
bubble by an aAbubble formed by two u/- propagator lines, for uj' ^ uj. In the 
following Fig. 12, we display some examples : 




Fig. 12 

(The corresponding "diagrammar" has been defined in Sect. 3 ; the summation, 

associated to each electron loop, is not displayed explicitly.) 
A characteristic feature of bosonization is that it reproduces only diagrams con- 
taining loops of at most two propagator lines. 

Introducing the effective interaction vertex 



UJ uj' UJ uj 1 UJ uj 1 UJ uj 1 




we can represent the diagrams reproduced as follows 




Fig. 13 
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The first class of diagrams leads to a contribution of order 1, the second one to 
a contribution of order y- . One finds that all leading order diagrams are repro- 
duced by bosonization, except for self-energy renormalizations of inner propagator 
lines. For short-range interactions, one can show by explicit calculations that the 
contribution of these self-energy renormalizations of inner propagator lines is zero 
- to the order considered in our calculation. For long range interactions, however, 
they could change the final result. 

One should remember that we use the effective action S n as an input of bosoniza- 
tion. Because of the linearization of the pieces of the Fermi surface contained 
in the boxes B w , this action describes the properties of the electron system cor- 
rectly at (unrescaled) momentum scales between ^ and . Therefore, eq . (p0|) 
reproduces the electron propagator only correctly for arguments |£ — rj\ of order 
between and ^ ( eq.(^) is written with respect to the rescaled system ). In 
this range, formula (|80[) holds to leading order in a (formal) expansion in powers 
off. 

To determine the electron propagator for larger arguments, we first have to cal- 
culate an effective (fermionic) action S m , m > n , describing the properties of 
the system at larger distance- and time-scales. Such an effective action S m can 
serve as a new input for the bosonization procedure. 



We now return to the calculation of the bosonic propagator of the interacting 
system. 

First, we consider systems with a short-range two-body interaction potential. 
For the non-interacting system, eq. (|74]) implies 



— — — In (k) . (81) 

kl + {v F ukf 

If one turns on interactions given by eq.(^), the propagator changes to become 



{27r) d+1 5 {d+1 \k+k') 



kl + (v Fn uky 



E 



\ A n 



-r n (k) 



(82) 



where 



g n (u ■ u') (kl + (v Fn ujk)(v Fn uk)^ + 
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') (~^0 + (vFn£k)(v Fn Uj'k) 



q n [—uj • oJ , , , . , . 

V /v u v " ' n k 2 + (u'k) 2 

For sufficiently small coupling constants, \g n (£u ■ uj') \ < g c <C 1 , the Neumann 

series converges, and one obtains 



(0 M (k)0 M (k')) n = (2-K) d+1 5^ d+1 \k + k') 



1 



kl + (v Fn Qk) 2 



1 + 



1 



jf(k;g r , 



(83) 

where f{k;g n ) is a bounded function in A;. From Sect. 2, we know that - for a 
system with short-range interactions - the coupling constants g n {uJ ■ u') tend to 
finite limits, as n — > oo. Thus, in the limit A n — > oo , the effects of interac- 
tions disappear, and the system is driven to the non-interacting (Landau liquid) 
fixed point with a propagator given by (|8~lD. For short-range interactions, the 



dependence on k± is irrelevant and suppressed by a factor 



Below, we shall see that, for sufficiently long-range transverse current-current 
interactions, the dependence on k± is significant ( i.e., singular, as \k±\ — > , 
with |fc | } \ k\\\ < \k±\\). It is not suppressed by an inverse power of the scale factor 
A n . We have to average over the variable k± in order to determine the effective 



dynamics of the bosonic degrees of freedom tp in the direction along 
Using eq.(PD|), we then calculate the propagator of the radial, quasi-electron com- 
ponents -0^1 . 

We start by studying two technically easier classes of systems with singular 
interactions. These systems correspond to examples (i)-(iii) of section 3 for which 
the polarization tensor U pcr has been calculated explicitly. 



i) Long-range, density-density interactions 



The set of coupling constants is given by eq. (|62]) , with g n V(\p\) = Tz=r\^S-\ a , a > 

K p P 

One obtains 

( 0[ u ](k) <p [u] (k>) ) n = (2n) d+1 5 {d+1 \k+k') 1 



kl + {v Fn uk) 2 
1 d (v Fn ujk) 2 



\ d - iy 2ix J k 2 + {v Fn uk) 2 



k 



g^v-Wf- l) + n?(fc) 
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where 11°° is defined in eq.([44|), with v F — > v Fn . 
Using that 

lim n°°(A:) — > cA 2 and lim Ii^(k) 



one verifies that 



(v Fn Ujkf 



kl + (v Fn iok) 2 



g^v-W-D + n^k) 



wk\2 
k I 



( 9n 


k 


\ ud-1 
t^p 


> 1 fc F A 


9n \ 


-11 k 


tup 


1 fcpA n 



+ c 



+ c' 



_fc_|2 

fco ' 



for | 



f or |& 



— 

fco 

fco 








For < a < 2 , the rhs is a bounded function of k. 

Under the assumption that all parameters V n (except the scale parameter A n ) 
characterizing the effective action S n tend to finite values, as n — > oo, it follows 
that, in the limit A n — > oo, the effect of interactions on the bose propagator 
disappears, as for the system with short-range interactions. 
Hence, in the limit A n — > oo, the electron propagator G tends to the standard 
LFL form. 



This result is due to the screening of the bare, long-range interaction po- 
tential V^dpl). In eq. (|84]) , the bare potential V^d^l) is replaced by an effective 
RPA— potential 



9n vf(k 0l k) 



9nV(\t\) 



1 + g n V(\±\W(k) gn'V-'m) + n°o(A;) 



In the static limit, 1^1 — ► 0, we obtain 

1 k 1 



lim g n vf(k ,k) 



:#t)- i (^) q + c 



i.e., in this limit, the effective potential is short-ranged. 
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ii) Tomographic Luttinger Liquid 



With eq . (|64j) , one obtains 

( <P{u](k) <P[u](k') )n = {2^) d+1 ^ d+1 \k + k') 



ko + (v Fn uky 



1 + 



9n_ 

2tt 



i -i 



As there is no dependence on the perpendicular momenta, one can immediately 
read off the effective dynamics for the modes ip^ propagating along [u] : 



2 Jnx[-^f,^f 



dkodku 



<PU(- k o,- k \\ 



1 + 



2tt 



(fco + {vFnk n ) 2 ) <£[L]Oo, fc|| 



The result for the quasi-electron propagator can be found with the help of 
formulae (A18) and (A19) of appendix A : 



where 



and 



e = (cei 



■ with 



The exponent 77 is given by 

7] 



(J d-l k e ifc ± a _ 



2f^) 2 



1 + 2 



'Sn\2 



This interacting system describes a Luttinger Liquid. 

The flow of the parameter Z n in eq.(p5[) can be derived by assuming that, along 
each direction u, the system is scale invariant : 

Recall that 

G(XnO ~ E e*'^' 4 G w ,„(0 ~ /dw e***-*" i- Gj , n (C) . (86) 



xV 
71 
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This equation and the assumption that the theory is scale invariant along each di- 
rection u implies the following matching-property for the Green functions G\ n (£) 

GijMQ « ^4,n+i(C) • (87) 

By inserting eq.([35"l) : 



(MCI 1 ^ M |C| 1+,? ' 

it follows that 

1 1 



J n+1 - rv 

~ or Z„ ~ 



Z n M* n (M n ) r i ' 

i.e., for this system, the residue Z n of the one-particle pole vanishes, as n — > oo. 



iii) Long-range, transverse current-current interactions 



9W 



Such interactions have been introduced in eq.(|66|), with g n V(\p\) = c n|^| a > c n 



3^ , and lead to the following bosonic propagator (for 



CO = CO 



'}) : 



1 - 



2 1 



Z-2 
A 



Xt 1 2n 2n i + 9n V(\*\) W n (k) k* + (v Fn uk) 2 



Co - 



k 2 



with 



d 



k 2 



(v Fn uky 



u 2 - 



{uokf 
k 2 



One finds that lim^oo n^(fc) = II* (/c) , with II* (k) defined in eq.( 



Using that lim^^H^k) ~ |-^| , one can see (cf. appendix C) that the contri- 
bution from the interactions to the propagator yields a singular /^-dependence, 
as \k±\ — > and (\k\\\, |Afa|) <C \k±\ . As discussed above, we then have to 
average over the perpendicular momenta kj_, in order to determine the (effective) 
action, , for the modes propagating along [to] . 
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The somewhat tedious calculations are deferred to appendix C. Here, we only 
describe our results. 

We obtain two different regimes, depending on the exponent a which characterizes 
the singularity of the interaction potential in momentum space by V'(ljpl) ~ : 

• For < a < d — 1 , the bosonization calculation yields 

G.AO ~ £V Mg( * 0+ ^ } Tl = c2 • (89) 



The function 8^ k ~)(£±) has been defined in eq.flS5|). The influence of the 



J,kp 

interaction on the electron propagator is suppressed by a factor A~^ d_1 ^ -Q! ^ 
as A n — > oo. 

From the matching condition (|8T|), it follows that 



Z n+ i ps Z n and Vp n +i ~ v Fn 
The system tends to a LFL. 



Kill . 



For d — 1 < a < 2 , we display the result in the region |£ | <S 



G u , n (g) ~ ^, fcF ; (^±) (^) e VF - Z n - 



exp ( -u l k F ck +a \K^zr\ 1+a ' ) , for |£||| > c„v Fn A 



q-(d-l) 



exp (-UnJ&.JM) , forc nWFn A° > > 



||||| ' iwi fe F ^ I 1 *!! I fe F 



for^i»|C||i»i 



(90) 

2 

The constant c n stands for 9 "If" , and m I5 m„ are positive constants depending 
on the dimension d of space and the exponent a (cf. appendix C). We have 
neglected all terms which are of lower order in an expansion in A" 1 than 
the leading terms displayed. Two interpretations are consistent with this 
result : 

1 . The parameters V n characterizing the effective action S n are all of or- 
der 1 - except for the scale parameter A n - , i.e., they do not flow 
under RG transformations. Then, the system tends to a LFL, as 



51 



A n — > oo : for finite arguments |£|||, the propagator G^n has the 
standard LFL form. It deviates from this form only for very large 

\d-l 

arguments ^> Fn k ™ — , where the results of a bosonization calcu- 
lation are not reliable (remember that, because of the linearization of 
the pieces of the Fermi surface contained in the boxes Q w , the results 
of the present calculation are only reliable for arguments |£||| smaller 
than A n ). 

2. The parameters of the set V n flow under RG transformations, i.e., they 
are functions of the scale parameter A n . 

The electron propagator deviates from the standard LFL form for finite 
arguments |£|||, if the condition 

C n V F nK ~ (j- ) 

is satisfied. Further, the requirement that the system is form- invariant 
under scale transformations implies the matching condition fl37|). This 
condition can be satisfied if 

d d—l — a. 

ck +a v Fn +a ~ const. . 

Then it follows that 

_ a(a — d+l) _ gd 

c n ~ An 1+a and v Fn ~ A n 1+a ■ 

The flow of the parameter Z n cannot be derived in the same way : for 
this, one would have to know subleading corrections to the argument 
of the exponential in eq. (|90|) . 

Hence the system displays "Non Landau Liquid" behaviour, as sug- 
gested in 0. 

The result of the bosonization calculation permits two consistent scenarios. 
In order to be able to decide which one is realized in the physical system one 
has to determine the flow of the parameters V n under RG transformations. 
Preliminary calculations indicate that the second case is realized (cf. also 

iBBB)- 
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For the critical value a c = d— 1 of the exponent a, we find for the electron 
propagator in the region |£ | <C ' 



G u , n (g) ~ ^, fcF ; (^±) (73—) e Z n 

Z7T 




- w /c F c n In 2 I ^|| I J , for I f || I > c n v Fn X r 



, for c n v Fn X n > |f || I > ± 



(91) 



where v is a constant of order one. 

As for a > d— 1, we can distinguish two possible scenarios, consistent with 
the results of the bosonization calculation : 

If the parameters c n ,VFn are not functions of the scale parameter A n , the 
propagator has the standard LFL form, as A n — > 00. 

In order to obtain a non-LFL behaviour, the product c n VF n X n of parameters 
must be of order one. Then, the matching condition (^71) requires that 
c n+ i ~ c n . By the first condition, it follows that Vp n ~ t-. 
Again, in order to decide which one of the two scenarios is realized, the flow 
of the parameters under RG transformations must be investigated. The 
calculations in |3], f|, |5], §j point in the second direction. 



Appendix A : Bosonization of the Electron Prop- 
agator in 1 + 1 Dimensions 

In this appendix, we review the bosonization of a relavistic 1 + 1 dimensional elec- 
tron system by using functional integrals. Our aim is to summarize the procedure 
for calculating the electron propagator. The general ideas of the bosonization 
technique have been presented in [[19], |(| . 



We want to investigate a 1 + 1 dimensional electron system in euclidean space- 
time, whose action is given by 

= Yd, tf(0 • (A.l) 
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Quantization is accomplished by using functional integrals. The field \& denotes 
a two-component Grassmann field, and \1/ := \l/*7 , where is an independent 
Grassmann field. Choosing the chiral representation of the 7-matrices, i.e. , 

7° = 01 7 1 = 02 7 5 = -«7°7 1 = 03 , 

the two components ^ ^ J of ^ are the antiholomorphic and holomorphic 

modes, resp., which are the euclidean analogues of left- and right-movers, resp. . 
We shall perturb the system by a current- current interaction of the form 

V($,V) = - 9 -Jd 2 U"(Of(0 , (A.2) 

where the one dimensional currents, j M , are defined by 

jt>($ = /:M7((W'M'U): ^ = 0,1 . 

This is the so called (massless) Thirrring model, or, in the context of solid state 
physics, the Tomonaga-Luttinger model. 

Before studying the interacting system, we establish the bosonization formulas 
for the free system. Identifying the fermionic current density j M with the bosonic 
expressions 

<-> f(<p;0 = -^-e^dMO , (A.3) 

'7T 



one can reproduce the Green functions for the currents 



J Zjb 

(A.4) 

where the action S(dtp) Q of the free, massless bose field p is given by 

S°(d<p) = Ifd^dMOdMO = \fMt)**MZ) ■ (A-5) 

Here, we have introduced the one form dip = d^ip d^ , and * denotes the 
Hodge star operation. We use standard notations of the calculus with differential 



^n this appendix, we write the action S((p) as a functional, S(d(p), of derivatives of the 
bosonic field. 
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forms. In |20| one finds a brief summary of the main definitions. ( Note that, 



in 



20 , we use slightly modified conventions; especially <p — > <p : — and 



In order to express the Green functions of the electric fields ( ^ ) in the 

V ^ J 

bosonic theory, one has to introduce disorder fields D(y, q) . 

For non-zero integers, q = {q^\ q( n >} , satisfying Y%=i Q = 0> we define the 

one form 

n 

n,,,(e) := E%), ? w(0 ( A - 6 ) 



1=1 



with 



n^ )j(?w (0 = nf* w (£)<r := -v^g (4) [*dA- 1 ^](0 

^7=? w 77 — W £ M^r = o, i 



Here, we use that 



d,d,A-\0 = -*«(£) ,i.e. , A^(0 = ~M£| • (A.7) 

Z7T 



In physics terminology, IL^ is the vector potential of a magnetic vortex of charge 
q at the point r] = (t]q,T]i) . The field strength is given by 

dU m (0 = ^q * 6{Z-ri) , 

If D is a two-dimensional domain then 



g , if G .D 
, else 



Locally on R 2 \ {rj} , the one form IL^ i? can be written as the derivative of a 
function 0£jj, g (£) , i.e. , 

U m (0 = da v , q (0 . (A.8) 

By identifying the euclidean space-time with the complex plane, £ = (£o)£i) l— ¥ 
£ = z£o + £i > the function a^ g (£) can be represented as 

= -^arg(£-r/) . (A.9) 
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Note that this function is not globally defined on C ; we identify the branch-cut 
with the negative, real axis. 



The expectation value of a disorder field D(rj, q) = rT?=i D(rj^\ gW) is given 

by 

I J r J ren 

On the r.h.s. of eq. ( |A.10| ), a multiplicative renormalization is necessary, in order 
to eliminate an infinite " self-energy" . 

Correlation functions involving disorder fields D(rj, q) and functional, J 7 (dip, if), 
depending on the bosonic field ip and its derivative dip are defined by 



D(r),q) T(dip,ip) ) B 



J Dip T(dip + Il M , ip + a M ) e 
J Vtp e- s (*e) 



-s(dip+n Viq ) 



ren 
(A.ll) 



(A.12) 



For the electron fields, the following identifications hold : 

(27r)3 

MO <— > - L r 

(2tt)3 
(2tt)3 

^(0 — 

(27r)3 

with the convention that, in a product composed of several if)' a s , we write the 
disorder fields to the left of all functionals depending on ip . 

The normal order of the exponentials is defined by 

. e i j d?z V (om . . = e ij <p$<p(om e |/^/(o(A+ m 2)" 1 /(?) ; (a. 13) 

where A = —d^d^ , / denotes a test function, and mo is a positive constant. 



£>(£,!) : 






£>(£,!) : 



















It is a staightforward calculation (cf. |19], P0fl ) to verify that 



n 

i=i i=i 
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Remark: We use the prescription that the function ( EKLi V^CCi) ) * s the 
analytic continuation of the function ( T IliLi V4- (&) ) > where T denotes the 
time-ordering operator for Grassmann fields, i.e., 



T 



^"(6) • ■ ■ I := sgn(a) ■ ■ ■ ^ a{n) ) , 



where a is the permutation such that > £°, 2 s > • ■ ■ > 
We shall make use of the following basic identities (cf. |19| ) : 







2 

nr=ic(^o)^ e _ ^ El ^" £l£jln( W ) , ifE"=i£j = , 

where the constant m can be choosen such that c(m ) = , 
1 



o , d 2 C nf, a (c)ng a (C) = 

Ea-^o+Ea- 1 ^) 



7T 

2 



L i=l 



i=l 



E M&-%l- E ln l6-0l- E ln l^-%l 

l<i<7<n l<i<j<n l<i<j<n 



and 



/ A nj? 2 (c) 5^(0 = o 



where the vortices located at the points £ := • • • , £ n ) carry the charge q — —\ 
and the ones located at the points rj :— (rji, ■ • ■ , rj n ) carry the charge 9=1. 

As an application, we calculate the Green function (ifj^iO^iv)) f° r the free 
system. Let £o > Vo- Then 



1 :( D(£, 1) : e^ (5) :: e"^^ : )° B 



f 2n 



^7T 



3 i-\/7r(a{,,( 1 f)-«e^(»?)) . 
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e -f[arg(j?-f)+arg(f-jj)] e -±ln|£-7)| . 



27T 



z e — : e * :: e 



1 



(A.14) 



2tt z(£ - 770) + (6 - T71) 
One observes that 

Similarly one obtains 

( m)Mv) ) = ^ —n ?T77 \ (A.15) 

and 

( ^1(0^2(7?) > = (v^(o^(»?)> = (v*(o#fa)> = ( mmiv) > = . 

These are the correct formulas for the free system ( compare, for example ap- 
pendix D of fll5], where the same conventions are used ). 



One can show that the same bosonization identities hold for the interacting 
system, with the interaction V given by eq.( [A.3|) (cf. p(J)- The on ly difference 



is that expectation values of bosonic operators are taken with respect to the 
interacting system, i.e. 

( ( . ) )V ■- f^t e -[ hf<% d M0dM0 + v(d v )] (. j ^ (A.16) 

J Zy 

where the functional V(ip) is obtained from the functional V{ip\ ip) by using the 
identities ( |A.3| ), i.e. , 

V(d<p) = £ J d 2 i dMOdMO ■ (A.17) 

We repeat the calculation of the Green function (-0| {C)^2 (v)) f° r the interacting 
system : 

(m)Mv)) v F = (1${<p;Z)M<p;v))b 



J ren 



, , ... IB 

'111 
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where 



Similarly, 



'2ix 



ren 



'2ix 



e -iarg(?-ry) e -±(l+f) In . 



Zy 



-ia,rg(£~ri) 



1 



2tt 



i +2 (?y 



(A.18) 



1 



e «arg(£-r?) 

2tt K-^| 1+/3 



(A.19) 



The exponent (3 describes the decay of the elecron propagator for large ar- 
guments. The dependence on the coupling constant is a characteristic feature of 
Luttinger Liquids. 



Appendix B : Consistency of bosonization formu- 
las 



In this appendix, we show on the example of ( Jr^i^) J\uAji) )° that the 
bosonization formulas (|70|) - (^) for the electron fields ip*, imply the bosoniza- 
tion formulas (|46| ) - (^) for the electron currents and j^. The generalization 
that this statement holds for arbitrary products of current densities evaluated in 
the non-interacting ground state is staightforward. 



Remark : 

The calculation of expectation values of products of current densities in the in- 
teracting ground state is organized as a perturbation expansion in powers of the 
set of coupling constants g n (uJ ■ 0') and the inverse scale parameter A" 1 . For sys- 
tems with interactions of the form given in eq. fl5"3|) (i.e. interactions which can 
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be expressed in terms of the current densities j^j and jr^i , cf. eq. ( p6| ) ) , such an 
expansion reduces the calculation to the task of evaluating products of current 
densities j^j and jr^i in the non-interacting ground-state - where our statement 
applies. 

By inserting eq. fl3~8|) , one obtains for the expression ( 3r wi ](£) Jh^Ajl) )° : 



e 



In order to simplify these expressions we can apply a " discretized" version of 
lemma (3.21) in jn| : 



Lemma : 

For a decomposition of the surface Sf~ l of the d-dimensional unit sphere (with 
d > 1) into congruent, quadratic patches with sides of length i we define the set 

of d-dimensional unit vectors cJj , % = 1, ...,N = Vol j^} 1 - pointing to the 
centers of these patches. Given a vector Bj G .M and some test function /(t^; £), 
one finds, in the limit A — > 00, the following asymptotic formula : 
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(B.l) 



where 



and 



zA Ald-l 
2 ' 2 1 



Here, A < Afcj? denotes an arbitrary momentum which can be chosen to be equal 
to kp. 



The proof of this lemma is analoguous to the proof of the continuous version of 
the lemma presented in appendix C of |]15| . 

By applying this lemma to the expressions appearing above and by inserting 
the bosonization formulas (|70|) - (|73|) for the electron fields ip^, we obtain the 
following formula : 

(^(Ojm^)) « -n^r- 1 ^ 2 Cd, ■ 

( dii^o^ii;^])^!!!^,^;^]) ) s0 ,j ) ( ^11(^0,^11; ^Lj^xiiC'm.^ii^ILj) ><?o ( ji ) ' 

(B.2) 

where Vu,||(£o' vL]) * s a short-hand notation for the bosonized expressions of 
the radial electron fields (£0, specified in eq. (|72|) . The explicit form of the 
radial electron propagators is given in eqs. (|A.14|) and (|A.15| ). 
It is a straightforward calculation to verify that eq. (B.2|) is identical to the fol- 
lowing formula: 

( 3[ U1 ](0 j[u 2 ](v) )° ~ 

= ' (R3) 
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where we use eq. (|A.7| ) for the propagator of the radial bose field ¥>L(£o, £||), an d 
eq.fl75|) to deduce the second part of the equation. 

Formula (|B.3|) reproduces the bosonization identity for the current density j u , cf. 
eqs. © -(H). 

Hence, we have verified our claim thus showing the consistency of the bosonization 
formulas for the electron fields and the current densities. 



Appendix C : Calculation of the Electron Propa- 
gator for a System with Long-Range, Transverse 
Current- Current Interactions 



In order to apply eq . (jSOD to the calculation of the propagator 

Gfc(£o - Vo, C|| -v\\) ■= - ( Vw(£o, £||) ^{Vo, v\\) ) l n 

of the radial electron field ip^ a , we have to determine the effective action S^(</>'L) 
of the bosonic field t^L propagating along the direction [a;]. 

The propagator ( <p[ u ](k) (p^k 1 ) ) n of the bosonized interacting system is 
given in eq. ([S8l) . It can be written as 

(%#)^#')>n = (2n) d +H^(k + k') 1 ^- Kfa(k) , (C.l) 

kl + {v Fn uk) 2 

with 









(1- 










(C.2) 



One observes that the function JC^(k) is bounded by < JCjL(fc) < 1 , for 
arbitrary k . 

In explicit calculations, expression (C.2) for JCf^Ak) is not convenient. We replace 
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it by 



l + g n V(\^-\)Ui(k) 



1 + 9»H\£\) 



W n (k) + 



k- 



2 1 (k F \d-l 



(C.3) 

which has the same asymptotics for small momenta k as )C^(k) . The function 
Unik) has been defined in eq. (|45|) , with vf — > VFn- 

The action S%((p[ u ]) of the bosonic field 0[ u ](k) which reproduces the propagator 
( |C.1| ) - but with the function )C(k) replaced by K.(k) - is given by 



We define 
with 



., , _ d a d+1 k0 H (-k) 

A JUxQ 



{kl + {v Fn ukf) (Kfa(k) 



-i 



(*£](*:)) 1 := 1+T^(k) 



(C.4) 
(C.5) 



2 1 (k F \d-l k o 



to 



2 _ (ukf 

1 



The effective action £j(<^L) °f t ne modes propagating along the direction [u] is 
deduced from the action S%((p[ u ]) by averaging over the momenta perpendicular 
to [oj] . One obtains : 



with 



- / dkn , 



dk {l <p [u] {-k) [(kl + {v Fn k\\f) (l+tf w] (feb,fc||))] [u] {k) 

(C.6) 



f#o,fc||) := -^J d^T^k) . 



kp J Q 

For small momenta (/so, fcn) and a large scale factor A n , the calculation yields the 
following asymptotic formula for t?^(k , k\\) : 



(d) *tt 



kp_ 1 

^+(^Fnfc||) 2 |fc | A^" 1 



fori > |c n A;o| > ^ 



a(d) 



fco + (fFnfc||) 2 



■r(^)r( 



V 1+a 



l7r ci-l2<J-2j + -L 



for-i > |c n A;o 



(C.7) 
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where d = 2, 3 is the dimension of space, a(d) is a positive constant of order unity, 
and the constants c n := -r£rtg n > and a, with < a < 2 , parametrize the 

k F 

interaction potential : g n V{\p\) = c n (j^) a ■ 

For 1 3> \c n k \ ^> , one observes that t n takes values between -j=x and 

\a-(d-l) 
n 

If a < d — 1, it follows that the interactions can be neglected in the limit A n — ► oo. 
If a > d — 1, the interactions become important as soon as fco is of order — or 
smaller. 

In the range ^> \c n k \ , eq.( |C.7| ) can be rewritten as 

k 2 



b(d,a)^r . * d , k F c^ a (^)^r^ fora>d-l 



2 7, 1 / A^, 1 g n |fco| \ 

K FC n in{ kp7r d-i 2 d-2 ) 
[ b(d,a)k F c n X a n ^ d -^ 



for a = d — 1 
for d - 1 > a > 



(C.8) 

For a < d — 1 , the effect of the interactions is suppressed by a factor d -i- a ■ For 
a > d — 1, one observes a singular ko - dependence of the function t. 



Using the methods presented in appendix A, the calculation of the radial elec- 
tron propagators is straightforward. Following eq. fl80|) , the large- distance asymp- 
totics of, e.g., the propagator G^ n is obtained by evaluating the expression 



Ci(C) Mo) > 



' fo;l 



(C.9) 

with C = (£o, ~~)- For each ray [ui], the vector potential and the function 



c.o 



are given in eqs.( |A.6j ) and ( |A.9j ) in appendix A. However, in contrast to 



appendix A, we here display the Fermi velocity vp n explicitly. Then, on the 
rhs. of eqs.( |A.6p and ( A. 9 ), there is an additional factor J— . Note also the 
supplementary factor ^v Fn in the exponent of the vertex operators. 

We again write the action S^X^L) as a functional, S^t^yiL), of the derivatives 
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One has that 



T/nM 



where the action S^(^jL) is defined in eg. ( |C.7| ). The contribution S^(H^q) of 
the disorder fields is given by 



and the contribution of the vertex-operators in eq.(C9) amounts to 



: : e 



exp 



TV I rf 2 k \p ik < -1] f VFn — - 

U 1 \k% + (v Fn k\\)* 1+tUk) 



(C. 



1) 



We display the behaviour of the fermionic Green function G^ n (Q in the following 
two regions : 

I = (C = (£0,-^4 : c„A^«|C|) , 
{ v Fn J 

II = { C = (£o, ■ ICI « c n \ a n and |£ | « 1-^1 1 • 



VFr, 



VFr, 



In region I , the propagator G^- n {Q is obtained by evaluating the integrals in eqs. 
( |C.10|) and ( |C.11| ) by using the asymptotic form ( |C.9| ) : 



( Ci(C) ^i(o) )1 



v 2tt 
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exp — Ujkp Cn +a A 



a-(d-l) 

1 + a \ 1 + a 

n 



ICI 1+a 



1+a 



a-(d-l) 



1 "Fr 



: COS 

||C| T+a <"' (,S 



^iarctan(|^|) 



1+a 



d arctan(|^^|) 



exp (— v k F c n In 2 \Q 



(i + Wl fe)ei 2 + (^neo 



e 2 + ( w ) 2 



kpc n 



2 > a > d- 1 



a = d — 1 



d - 1 > a > 
(C.12) 



where 



it, 



c(d, a) 



7T 



1 + a sinTr^; 



r( 



d 



and 



c(d, a — d — 1) - 
a 



Remark : As it stands, the limit of the propagator for a { d — 1 is singular 
and does not reproduce the formula for a = d — 1 . The reason is that, for 
a > d— 1 , there is a subleading term, not displayed in eq. ( |C.12|) , which becomes 
leading for a = d — 1 . Idem for a] d — 1 . 



In region II, one finds 



for d - 1 < a < 2 

(Ci(C)^i(o) >|i 



. —% 

•2tt 



)e 



sarg(2£ + — ^) 



- w nra=r 



Kill 



forc ri ^ n A«»|e|||»^^ i , 
for " F " A " 



hp 



(C.13) 
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• for < a < d - 1 



Ci(0 ^i(o) >i 



- l7Tr) e u Fn —— - , f or |£|||»— . 



•2vr' 



(C.14) 



Remark : In eqs. flC.13 ) and ( C.14 ), we have omitted sub leading terms in an 



expansion m y- 



An 
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